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Abstract 

The consistency of loop regularization (LORE) method is explored in multiloop calculations. A 
key concept of the LORE method is the introduction of irreducible loop integrals (ILIs) which are 
evaluated from the Feynman diagrams by adopting the Feynman parametrization and ultraviolet- 
divergence-preserving(UVDP) parametrization. It is then inevitable for the ILIs to encounter the 
divergences in the UVDP parameter space due to the generic overlapping divergences in the 4- 
dimensional momentum space. By computing the so-called a/37 integrals arising from two loop 
Feynman diagrams, we show how to deal with the divergences in the parameter space with the 
LORE method. By identifying the divergences in the UVDP parameter space to those in the 
i ' subdiagrams, we arrive at the Bjorken-Drell's analogy between Feynman diagrams and electrical 

circuits. The UVDP parameters are shown to correspond to the conductance or resistance in the 

J> 1 electrical circuits, and the divergence in Feynman diagrams is ascribed to the infinite conductance 

or zero resistance. In particular, the sets of conditions required to eliminate the overlapping 
momentum integrals for obtaining the ILIs are found to be associated with the conservations of 
electric voltages, and the momentum conservations correspond to the conservations of electrical 
currents, which are known as the Kirchhoff 's laws in the electrical circuits analogy. As a practical 
application, we carry out a detailed calculation for one-loop and two-loop Feynman diagrams in 
the massive scalar <fr theory, which enables us to obtain the well-known logarithmic running of 

j> . the coupling constant and the consistent power-law running of the scalar mass at two loop level. 

Especially, we present an explicit demonstration on the general procedure of applying the LORE 
method to the multiloop calculations of Feynman diagrams when merging with the advantage of 
Bjorken-Drell's circuit analogy. 
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I. INTRODUCTION 



Quantum Field Theory(QFT) is the most successful theory for understanding the micro- 
scopic world in elementary particle physics, nuclear physics and condensed matter physics. 
However, when carrying out any calculation beyond tree level in the framework of perturba- 
tion treatment of QFT, one would encounter the infinities in Feynman integrals, coming from 
contribution from large momenta and are usually called ultraviolet (UV) divergences. Thus, 
the QFT becomes well-defined only when it can be regularized and renormalized properly. 

Nevertheless, the widely-used regularization methods are known to have some limitations. 
For instance, the Pauli-Vallars regularization method is very useful in the calculation of 
quantum electrodynamics(QED), but it fails in non-Abelian gauge theory as it explicitly 
destroys the non-Abelian gauge invariance. The dimensional regularization can preserve 
the gauge symmetry explicitly and has been useful in the computations for gauge theories, 
such as QED and QCD of the standard modelpj. Despite its great success, it has been 
known [j], Q that the spinor matrix 75 and chirality cannot in principle be well defined in the 
extended dimensions. Also it has trouble in applying directly to supersymmetric theories 
which depend dimension of space-time, and moreover it cannot keep track of the divergence 
behaviors (quadratic and above) of original integrals in the Feynman diagrams. So it is not 
useful in some calculations in effective field theories and chiral dynamics where we need to 
isolating the quadratic divergences for understanding the dynamical symmetry breaking and 
restoration. 

Thus it is desirable to develop an alternative new regularization scheme which possesses 
the basic properties: being well defined in 4-dimensional space-time, preserving the gauge 
symmetry and Lorentz symmetry, keeping track of the divergent behaviors of original theo- 
ries, making the pratical calculations as simple as possible and applicable to both underlying 
and effective QFTs as well as supersymmetric and chiral QFTs. 

Recently, a new regularization method proposed by one of us[l, 0] can satisfy all of 
the properties mentioned above and checked carefully with explicit calculations for many 
applications at one-loop level. For convenience, such new regularization is called the Loo 
Regularization since its prescription acts on the so-called irreducible loop integrals(ILIs) 
5J. For short, here we may use 'LORE' as an abbreviation of the loop regularization. It 
has been proved with explicit calculations at one loop level that the LORE method can 
preserve non-Abelian gauge symmetry [|3| and supersymmetry 0] . It can provide a consistent 
calculation for the chiral anomalyj3] and the radiatively induced Lorentz and CPT-violating 
Chern-Simons term in QEDQ as well as the QED trace anomaly|9]. It allows us to derive 
the dynamically generated spontaneous chiral symmetry breaking of the low energy QCD 
for understanding the origin of dynamical quark masses and the mass spectra of light scalar 
and pseudoscalar mesons in a chiral effective field theory (Io|, and also to investigate the 



chiral symmetry restoration in a chiral thermodynamic model[ll[|. In particular, it enables 
us to consistently carry out the quantum gravitational contributions to gauge theories with 
asymptotic free power-law running 12Ml4 . 

It has been analyzed in ref . 0] that the LORE method can straightforwardly be generalized 
to higher loop calculations with an explicit demonstration on the general two loop integrals, 
i.e., the so-called a/?7 integrals. In fact, our general proof for the consistency of loop regu- 
larization via the a(3j integrals was just following the same procedure which was adopted by 



't Hooft and Veltman|l] to demonstrate the consistency of dimensional regularization. Since 



the LORE method has been realized in four dimensional space-time without modifying the 
original theory, its consistency cannot be proved in the Lagrangian formalism to all orders, 
thus it is useful to develop a diagrammatic approach to make such a general proof. For that, 
we shall make explicit multiloop calculations to show its consistency, figuring out a gen- 
eral procedure for practical calculations, which is a further motivation in our present work. 
We are going to show in present paper that the evaluation of the irreducible loop integrals 
(ILIs) from Feynman integrals by adopting the ultraviolet divergence-preserving (UVDP) 
parametrization naturally leads to the Bjorken-Drell's circuit analogy between Feynman 
diagrams and electric circuits. As a consequence, when merging the LORE method with 
the Bjorken-Drell's circuit analogy, we arrive at the interesting observation that there is 
the one-to-one correspondence between the divergences of the UVDP parameters and the 
subdiagrams of Feynman diagrams, which enables us to extend the procedure to higher loop 
Feynman diagrams in a more general and systematic way. 

The key concept in the LORE method is the introduction of the ILIs which are ob- 
tained from the Feynman diagrams by using the Feynman parametrization and the UVDP 
parametrization. A crucial point in the LORE method is the presence of two energy scales. 
They are introduced via the string-mode regulators in the regularization prescription acting 
on the ILIs. It has been shown that the two energy scales play the roles of the ultraviolet 
(UV) cut-off and infrared (IR) cut-off to avoid infinities without spoiling symmetries in the 
original theory^, It is then inevitable to encounter the UV divergence in the UVDP 
parameter space due to the generic overlapping divergences. To be more explicit, we carry 
out a calculation for the general integrals, the so-called a/?7 integrals, arising from two loop 
Feynman diagrams, and show how to deal with the divergences in the UVDP parameter 
space by applying the LORE method. By identifying the divergences in the UVDP parame- 
ter space with those in the subdiagrams, we naturally arrive at the Bjorken-Drell's analogy 
between Feynman diagrams and electric circuits, where the UVDP parameters are found to 
be associated with the conductance or resistance in electric circuits. A detailed description 
on circuit analogy is given in the book by Biorken and Drell 15] in which the circuit analogy 
was originally inroduced to study the analyticity properties of Feynman diagrams from the 
causality requirement. In our present paper, we observe that the sets of conditions required 
to eliminate the overlapping divergent momentum integrals for evaluating the ILIs is analo- 
gous to the conservations of electric voltages in the loop, and the momentum conservations 
to the conservations of electric currents at each vertex. These equations are known as the 
Kirchhoff's laws in electric circuits. In particular, it is noticed that the divergence in Feyn- 
man diagrams corresponds to an infinite conductance or zero resistance in electric circuits. 
By adopting such an analogy, we perform a detailed calculation for one- and two-loop Feyn- 
man diagrams in the massive scalar 4 theory, and meanwhile we explicitly demonstrate 
the general procedure for applying the LORE method to multiloop calculations of Feynman 
diagrams. 

We would like to emphasize that our motivation is not just for figuring out a much simpler 
regularization scheme, but for finding out whether there exists in principle a regularization 
scheme which can overcome some shortages and limitations in the widely-used regularization 
schemes. Meanwhile, we expect that such a regularization scheme must also be practical and 
as simple as possible. In fact, for the one loop calculation, the LORE method is really simple. 
For the higher-loop calculations, the procedure and calculation in the LORE method are 
not as concise as the ones in the dimensional regularization, but our treatment here makes 
the overlapping divergent structure as well as its divergent behavior physically manifest. 
To be more precise, the divergence structure for a diagram includes the overall quadratic 
or logarithmic divergence and the divergences in the subdiagrams, as well as corresponding 



subtraction diagrams. Actually, the simplicity of the dimensional regularization is at the cost 
of three essential limitations: (i) the definition of 75 in theories beyond 4 dimension, (ii) the 
requirement of exact dimension of original theories, (iii) preservation of quadratic divergence 
in original theories. To overcome such limitations in the dimensional regularization is our 
main purpose to look for a possible alternative consistent regularization scheme, which will 
be helpful for understanding deeply the applicability and consistency of QFTs. In this sense, 
the LORE method has been a step forward, as already shown in 0-Ei at one loop level. 
It is worthwhile to go further and make an explicit check at higher loop level, which is our 
main goal in the present paper. Of course, for QFTs without three limitations in principle 
mentioned, the dimensional regularization scheme remains a powerful and simple one for a 
practical calculation. 

The paper is organized as follows. In Sec. II, we briefly outline the LORE method and the 
concept of ILIs at one loop level. In Sec. Ill, a particular contribution of two-loop vacuum 
polarization diagram in QED is examined and show the general structure of overlapping 
divergences. It is then unavoidable to encounter the UV divergences hidden in the UVDP 
parameter space. In Sec. IV, we apply the LORE method to the general a/?7 integrals with 
a = 7 = l,/3 = 2, and explicitly show how the LORE method can appropriately regularize 
the UV divergences either from the original loop momenta or from the UVDP parameters. 
In Sec. V, we show how the evaluation of ILIs and UVDP parametrization naturally merges 
with the Bjorken-Drell's electric circuits analogy. In Sec. VI, The Bjorken-Drell's electric 
circuit analogy of Feynman diagrams allows us to analyze the origin of UV divergences 
contained in the UVDP parameter space, and to figure out the one-to-one correspondence 
of divergences between subdiagrams and UVDP parameters. In particular, the divergences 
in Feynman diagrams is shown to correspond to infinite conductances or zero resistances in 
electric circuits analogy. In Sec. VII, the LORE method combining with the Bjorken-Drell's 
analogy shows the advantage in analyzing a complicated overlapping divergence structure 
of Feynman diagrams. As an explicit illustration, the case with a = /3 = 7 = lof the 
general a/37 integral is discussed in detail and all the harmful divergences cancel exactly. 
As a practical application of all the machinery privously introduced, we carry out in Sec. 
VIII a detailed calculation of two loop contributions in the massive scalar 4 theory. Some 
additional quadratic corrections to the scalar mass are obtained, and leads to a power-law 
running. Based on the general analysis and explicit calculations, we arrive at in Sec. IX the 
general procedure of applying the LORE method to high-loop calculations. Our conclusions 
and remarks are presented in Sec. X. 



II. CONCEPT OF ILIS AND BRIEF OUTLINE ON THE LORE METHOD 

We start from the fact that all Feynamn integrals from the one-particle irreducible (1PI) 
graphs in 1-loop can be written, by using Feynman parametrization, in terms of the following 
sets of loop integrals, 

f d 4 k 1 

2a ~ J (2^ (k*-M 2 ) 2 +« ' a--l,0,l,2,--- , 

for scalar type integrals and 

I" d k k^ky 

f d k k^k v k p k a _ _i n 1 9 (\\ 



for tensor type integrals, where the number (—2a) in the subscript labels the dimension 
of power counting of energy momentum in the integrals. Thus two special cases a = — 1 
and a = correspond to the quadratic divergent integrals (J 2 , h/iu-) an d the logarithmic 
divergent integrals (J , io/^-)- Note that the mass factor M 2 is in general a function of the 
Feynman parameters and the external momenta Pi, namely M 2 = Ai 2 (m 2 ,pl, ■ ■ ■ ). 

The above loop integrals are the so-called one-fold irreducible loop integrals (ILIs)0], 
which can be generalized to the n-fold ILIs evaluated from n-loop overlapping Feynman 
integrals of loop momenta k%(i = l,2,---n). In general, the n-fold ILIs are defined as the 
loop integrals in which the overlapping momentum factor (hi — kj + pij) 2 (i ^ j) originally 
appearing in the overlapping Feynman integrals has already been eliminated. It has been 
shown that any loop integrals can be evaluated into the correspondinglLIs by using both 
the Feynman parametrization and the UVDP parametrization methods[3|. Note that in the 
procedure of evaluating the ILIs, the algebraic computing for multi-7 matrices involving 
loop momentum k 1 such as ^7 M ^/ should be carried out first and expressed in terms of the 
independent components: 7^, a^ u , 757^, 75. 

The concept of ILIs is crucial in the LORE method. To see that, let us briefly examine 
the vacuum polarization in the non-abelian gauge theory. We begin with the following 
lagrangian in gauge, 

C = Mi^D, - m)^ n - - A F; V FH V - hdPAff + PifcDtf; (2) 

with 

= 9/j,A^ — d^A^ — gfabcA^AI 
D^ n = {d^ + igT a Al)i; n , (3) 

where £ is a gauge parameter. ip n , A^ and 77 are fermions, gauge bosons and ghost fields, 
respectively. T a are the generators of gauge group and f a b c the structure constants of the 
gauge group with [T a , T b ] = if a b c T c . The vacuum polarization corresponds to the self- 
energy diagrams of gauge boson, which contains the quadratically divergent integrals, the 
most divergent behavior in all of the Green functions in one-loop. Here we give the final 
results carried out by using the usual Feynman rules in the general £ gauge. The details 
of the calculation can be found in ref. |3j. The explicit expressions for the gauge boson 
self-energy diagrams are given, in terms of the ILIs, as follows: 

^i f J ab = -g 2 4N f C 2 6 ab [ dx [ (2a 2 - l)h{m) 9liV 

Jo 

+2x(l - x)(p 2 g^ u - ptf>v)Io{m) ), 
for the fermion loop contribution to the gauge self-energy diagram, and 

n^ )afe = g 2 C x b ah {p 2 gilv - p^p v ) [ dx{[l+ 4x(l - x)\ J 

Jo 

+ [ ( l + 6x(l -x)(a + 2) -3a ) J - 2x(l - x) ( 1 + 12x{l - x) )p 2 J_ 2 ] 
3 

+-A 2 a_ 2 x(l - x) p 2 I-2 } 

+9 2 C 1 5 ab / dx { 2( 2a 2 - 1 )I 2 + A(a - 1) p^Pv x{\ - x) p 2 J_ 2 }, (4) 
Jo 



for the gauge boson and ghost loop contributions to the gauge self-energy diagram, where p 
is the momentum of the external gauge boson, Nf the number of fermions flavors, A = 1 — £, 
facdfbcd = Ci^ab and tr (T a T b ) = C 2 5 &. We have also used the following definitions from 
the general Lorentz decomposition 

1 1 

hfiv — 0-2 I2 9fiu, Ionv — ^ a -^0 9iivi I-2nv = ^ a -2 1-2 9\j,vi (5) 

where i_2 and 1 - 2vLV are convergent integrals with a_2 = 2/3. Note that IT^) ab depends on 
the gauge parameter £. This is because the Green's functions are gauge dependent while 
only the S- matrix elements are gauge independent. However, current conservation implies 
that Il ( i!J ab and njfj° 6 have to satisfy the Ward Identities p»Il$ ab = and pTlJfJ ab = 0. 
Notice that the first line of IllIJ ab and in the last line of IL^; contain both quadratically 
and logarithmically divergent integrals which might violate gauge invariance. Only with the 
following consistency conditions 

a 2 = 1/2, a = 1 (6) 

then the gauge invariance can be preserved. 

Nevertheless, from the naive analysis of Lorentz decomposition and tensor manipulation, 
one gets by multiplying g^ v on both sides of Eq.(j5]), 

9^l2^ = h + M 2 Io = Aa 2 I 2 , i.e. a 2 = l/4 + M 2 I /h, 
g^Ionv = I + M 2 I- 2 = a J , i.e. a = 1 + M 2 I- 2 /I , 

which leads, without using any regularization schemes, to the following relations 

1 1 

h/jv = 7^ h + -^g^uM 2 1 , 

1 1 i 

hnv = -^g^u Io + g^uM 2 I- 2 = -g^ u I - 32^2 0*"" ( 7 ) 

Clearly, the above naive relations for the divergent ILIs will destroy the gauge invariance. 
The reason is that in the divergent integrals which are generally not mathematically well 
defined without using proper regularization scheme, the tensor manipulation and integration 
do not commute with each other, so the result for divergent integration is not consistent 
in general. Thus in order to obtain a consistent result, one has to adopt a regularization 
scheme to make the divergent integrals well-defined. To see this, consider the time-time 
component on both sides of the relation for the quadratic divergent ILIs in Eq.([5]) 
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The Wick rotation will turn the four-dimensional energy momentum into Euclidean space 
and integrating over the zero component of energy momentum on both sides, we get 

d 4 k 1 r d 3 k r „ 1 
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d 3 k 1 
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~i I 7^77 2 - 1 tan' 1 ( k Q /y/\s? + M 2 
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for the right-hand side, and 



<2 oo 



d k Jvq , f d k f j . 



dk - 



(2tt) 4 (k 2 + M 2 ) 2 J (2vr) 4 J u (^ + k 2 + M 2 ) 2 

_■ / d3fc [jl f 1 k 2 + A4 2 
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2 J (2tt) 3 2Vk 2 + M 2 2 

for the left-hand side. Note that the above integration over fc Q is convergent, and should 
be safe for any algebraic manipulation. When comparing the results with both left and 
right hand sides, we obtain 02 = 1/2 which agrees with the consistency condition for gauge 
invariance. We then come to the conclusion that the general relation between the tensor-type 
and scalar- type quadratically divergent ILIs with 02 = 1/2 must be the exact consistency 
condition. 

We would like to emphasize that the above demonstration for obtaining the consistency 
condition 02 = 1/2 between the quadratically divergent ILIs has nothing to do with any 
regularization schemes. Nevertheless, the drawback here is that it is obtained only for one 
of the Lorentz components rather than for the whole covariant Lorentz tensor. Thus it 
is necessary to look for a proper regularization scheme which can realize the consistency 
condition in a covariant way with the well-defined divergent integrals. Meanwhile, it should 
also preserve the original divergent behavior for both quadratical and logarithmic divergent 
integrals. Actually, it has explicitly been proved ^ that the LORE method does lead to 
the consistency conditions with 02 = 1/2 and ao = 1. A simple regularization prescription 
operating on the ILIs has been realized in four dimensional spacetime to satisfy the criteria 
mentioned in the introduction. 

The regularization prescription of the LORE method is as follows: Firstly rotating the 
momentum to the four dimensional Euclidean space, then replacing the loop integrating vari- 
able k 2 and the loop integrating measure J d A k of the ILIs by the corresponding regularized 
ones [k 2 ]i and J[d A k]i: 



k 2 -> [k 2 ^ = k 2 + M l 



1 > 



N 

[ d 4 k^ I ld 4 k} t = lirn J^cf / d% (11) 

J J N,M l l=0 J 

where Mf (I = 0,1, • • • ) may be regarded as the regulator masses for the ILIs. The 
coefficients cf and the regulator masses are chosen to satisfy the following conditions: 



N 

lim^cf \M 2 r = (n = 0,l ) ...), (12) 

where the notation \im NM 2 denotes the limiting case lim A r M 2^ 00 . The initial conditions 
Mq = and Cq = 1 are taken to recover the original integrals in the limit Mf — > 00 



(1 = 1,2,... )• 



With the above regularization prescription, we have shown that the regularized 1-fold 
ILIs satisfy the following consistency conditions jijj]: 

R 1 R R 1 R 

^2iiv = 7,9p,v T-2 , Iluvpo = o(9(J.v9pa- + 9pp9va + 9pa9pv) ^2 



^Opu ~ a 9\iv Iq, ^Opupa ~ o/l ^9pv9pa + 9pp9ua + 9po9pv) Iq ■ (13) 



l 2p,v ~ 2 y ^ v 2 j ± 2fiupc7 ~ gKyp.vypv ~ yppyua ^ ypaypvj * 2 > 

I tR tR _ J_ 

^y/if 1 o ' 1 o/nypa 24 ' 

which are actually the necessary and sufficient conditions to preserve the gauge symmetry 
in QFTs. Here the superscript "R" denotes the regularized ILIs. Note that the dimensional 
regularization scheme also leads to a 2 = 1/2 for M 2 ^ and a$ = 1, while the resulting 1 2 is 
suppressed to be a logarithmic divergence multiplying by the mass scale M 2 , thus it goes to 
vanish 1^ = when Ai 2 = 0. This is the well-known fact that the dimensional regularization 
does not preserve the quadratic divergent behavior of the original loop integrals. 
As the simplest solution of Eq. ( TT2|) . take the string- mode regulators 

Mf = + IM 2 R , (14) 

with I — 1, 2, • • • , then the coefficients cf are completely determined to be 



-O'tat 77177' (15) 



' v ~' (N-i)w: 

where Mr may be regarded as a basic mass scale of loop regulator. When applying the 
above prescription and solution to the ILIs, the regularized ILIs in the Euclidean space-time 
are generally expressed as follows: 

N 



1-2* = i( 



-l) a lim V cf [ 



d^k 



., (27r) 4 (k 2 + M 2 + M 2 ) 2H 
tR _ j(_i)a H rj / d * k k^k u k p k a 

i-2*»»pa ~ H ij y (27r) 4 (A; 2 + M 2 + Mf) A+a { ' 

For the regularized quadratically and logarithmically divergent ILIs I-f and J(f , we have 
shown that they have the following explicit experssionsj3|: 



$ = 7tA{M c 2 - /*> ^ - lw + 1 + y 2 (4~ 2 )]} 



16tt 2L c ^ /x 2 ' ' " ZV M 2 



i M 2 11 2 
$ = T^[ln-f-7. + yo(^)] (17) 
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with [i 2 = /j 2 + M 2 , and 

AT AT 

7™ = lim{ ^ cf In / + ln[ ^ cf / In / ]} = 7s = 0.5772 • • • , 

z=i i=l 

y (x) = / da , y^x) = 

Jo x 

V2{x) = Vo{x)-y 1 {x) : limy^rr) 0, i = 0,1,2 (18) 

A 7 

M 2 ee lim M 2 R \^cf(l In/) = lim Ml/lniV 

/— 1 



which indicates that the fi s sets an IR 'cutoff' at M 2 = and M c provides an UV 'cutoff'. 
For renormalizable QFTs, M c can be taken to be infinity (M c — > oo). In a theory without 
infrared divergence, /i s can safely be taken as fi s = 0. In fact, by taking M c — > oo and 
fi s = 0, we recover the initial integral of the theory. Also by taking Mr and N to infinity, we 
arrive at a regularized theory which becomes independent of the regularization prescription. 
Note that the function yo{x) with x = fi 2 /M 2 is actually the incomplete gamma function, 
which has the property: yo{x) — » at x — > (i.e., in the limit M c — > oo). In comparison 
with the dimensional regularization, there is a correspondence: In ^f- — > ~ with M c — > oo 
and e — > 0, which indicates that the function yo(x) approaches to zero much faster than 
the polynomial of e in the dimensional regularization. This can be seen explicitly from the 
expression: yo(x) ~ x ~ — >■ in the limit M c — > oo and £ — > 0. 

We would like to point out that the prescription in the LORE method looks very similar 
to the Pauli-Villars prescription. Nevertheless, the basic concept is quite different as the 
prescription in the LORE method is acting on the ILIs rather than on the propagators in 
the Pauli-Villars scheme. This is why the LORE method can preserve non-Abelian gauge 
symmetry, while the Pauli-Villars regularization can not. In this sense, we would like to 
emphasize that the concept of ILIs is a crucial point in the LORE method to realize the 
interesting symmetry-preserving regularization scheme. In particular, the introduction of 
two intrinsic energy scales without spoiling symmetries of original theory is an advantage 
in the LORE method to avoid the infinities of divergent Feynman integrals [3- 14]. For 



the effective theories, the intrinsic UV 'cutoff' scale M c plays the role as the characteristic 
energy scale below which the physics can be well described by the effective quantum field 



theory [10H14 



III. OVERLAPPING DIVERGENCES AND UVDP PARAMETRIZATION 

It is well-known that for Feynman diagrams beyond one-loop order, a new feature, over- 
lapping divergences. It happens when two divergent loops share a common propagator. To 
illustrate this, consider one particular contribution to the photon vacuum polarization at 
two-loop order of quantum electrodynamics (QED) (see Fig. 1) 




FIG. 1: 

Here we may follow the argument in the textbook [16(]. As discussed in the usual texts 
of QFTs, the UV divergences in Fig. [TJ can arise from several regions of momentum spaces. 
One divergent contribution comes from the region where there is a large momentum passing 
through the left subdiagram. This means that the three points x, y, and z in position space 
are very close together, while the point w is farther away. In this region one can think 
that the virtual photon gives the corrections to the vertex x. Plug the divergent part of the 



one-loop vertex corrections into the rest of diagram and integrate over k\. We get expression 
identical to the one-loop photon vacuum polarization correction multiplied by the additional 
logarithmic divergence, as shown in Fig. [2j Obviously, a similar divergent contribution to 
the diagram in FigJT] arises from the region with a large momentum passing through the 
right subdiagram as shown in Fig. [2j It is manifest that the log 2 A 2 term comes from the 




~ a(g fJlU p 2 — p^)(log A 2 — logp 2 ) ■ alogA 2 
FIG. 2: 

region where both ki and k 2 are large. While the log p 2 log A 2 term results from the region 
where k 2 is large but k\ is small, another such a term arises from the region where ki is 
large but k 2 is small. The terms like log p 2 log A 2 are called nonlocal or harmful divergences 
as such terms cannot be canceled by the ordinary substraction scheme by introducing the 
corresponding two loop counterterms in the Lagrangian. 

It is then expected that these harmful divergences are canceled by two types of countert- 
erm diagrams. First, we can build diagrams of order a 2 by inserting the order-a countert- 
erm vertex into the one-loop vacuum polarization diagram (see Fig. [3]). Such two diagrams 




FIG. 3: 

should cancel the harmful divergences as shown in Figj2] Once these counterterm diagrams 
are added, the only divergence left is exactly local and can be canceled by the two-loop 
overall counterterm, which is diagrammatically represented in Fig. SJ 




FIG. 4: local counter term 

The discussion given above is a general description in the ordinary textbooks and there is 
no problem in principle. However, in the practical calculation, one actually meets some con- 
ceptional problems. There is no doubt that we have to integrate over two loop momentums 
k\ and k 2 one by one. Suppose that we first integrate over the loop momentum k\, which 



means that we integrate over the left subdiagram corresponding to the left vertex insertion. 
Then we integrate over the loop momentum k 2 , which corresponds to the overall divergence 
of the whole diagram as its divergent behavior is easily found to be quadratic from the 
simple power counting. But the problem arises when we look into which loop momentum 
integral represents the right subdiagram corresponding correction to the right vertex. This is 
because we have already integrated over both loop momenta in the diagram with the above 
procedure. It appears that we have nothing to do with it. Actually, when carrying out the 
calculations by using the Feynman parametrization and UVDP parametrization to combine 
the momenta in the denominator, we will find that, besides of the divergences coming from 
the integral of the two loop momenta k\ and k 2 , the integrations over the UVDP parameters 
are also logarithmically divergent, which exactly reproduce the divergence behavior of the 
vertex correction at one-loop order. This observation makes it clear that the integration 
of right subdiagram is "hidden" in or transformed into the parameter space with the usual 
procedure of dealing with the two-loop overlapping diagrams. 

Thus the next immediate question is whether, given a divergence in the UVDP parameter 
space, we can find out the origin of this divergence in the original Feynman diagrams. Our 
answer is positive. This is actually the main purpose in our present paper. We shall show 
that there is an exact correspondence between the UVDP parameter integrals and those 
from the original loop momenta. The key conceptual tool for arriving at this conclusion is 
the observation of the Bjorken-Drell's analogy between the Feynman diagrams and electrical 
circuits, which will be demonstrated below. 

Before proceeding, it is interesting to note that all the overlapping divergent integrals (in- 
cluding scalar-type and tensor-type) of two-loop Feynman diagrams in QED can be reduced 
to the following two types of integrals by adopting the Feynman parametrization: 

/d^ki f d 4 k 2 1 
(2^ J (2vr) 4 Wl - - ml)[(fci - h + p) 2 - m§] ' (19) 



d^h r d^k 2 

' ,2> ~ 1 (2^ J (2vr) 4 {k 2 - m\){k 2 - m|) 2 pi - k 2 + p) 2 - mj ' 



(20) 



where mf are in general the functions of the external momenta p and Feynman parameters. 
Such integrals are actually the two special cases of the general a/?7 integals[l[. Therefore, it 
is useful to make a general discussion and analysis on the regularization and renormalization 
for the general aftj integrals. 

In order to avoid the complication involving the reducible loop integrals and tensor-type 
integrals, we may consider only scalar-type ILIs. As pointed out in Ref. [jj by 't Hooft and 
Veltman, a general two-loop order Feynman diagram can be reduced to the general a/37 
integrals of the form: 

d^ki f d 4 k 2 1 ^ 



aM J (2tt) 4 J (27ry(k 2 1 -m 2 l ) a (k 2 -m 2 2 y[(-k l -k 2 + P y-mlp' 

We shall focus on the problem how to disentangle the overlapping divergences with the 
LORE method. Especially, we will show how to deal with the divergences contained in the 
UVDP parameter space caused by the overlapping structure. With the advantage of the 
UVDP parametrization in evaluating the ILIs from Feynman diagrams, we naturally arrive 
at the Bjorken-Drell's analogy between Feynman diagrams and electrical circuit diagrams. 
This powerful tool gives us exact one-to-one correspondence of the divergences between 
the parameter space and the subdiagrams, which can explicitly be demonstrated. As a 
consequence, it straightforwardly leads to the important theorem on the cancelation of 



harmful divergences. To realize those goals, it is enough to keep track of only the overlapping 
divergences, such as the terms ■ log -^p and log ■ log 3%. For the harmless divergences 
and finite terms, they can be either absorbed into the two-loop overall counterterms or kept 
in the final expression. 

From the general form of Eq. (1211) . one can easily recognize that there are in general 
one overall integral a(3^ and three subintegrals (a (3, f3 r y and 7a), represented diagram- 
matically as the following three corresponding subdiagrams (See Fig.Qand Fig.(j6])). The 
corresponding counterterm diagrams are shown in Fig.^7]), which are generally needed for 
the cancelation of the harmful divergences. 



&4, a, v\ 




FIG. 5: af3j diagram 




FIG. 6: subdivergences 




FIG. 7: counterterm diagrams 



By power counting, it is easy to see from Eq. ([2"Tj) that there are only two cases which 
involve the overlapping divergences, i.e., (1) a + /3 + 7 = 4, and (2) a + f3 + r y = 3. Other cases 



with a + (3 + 7 > 4 contain only harmless divergences and the overall integral is convergent. 
Thus we shall only discuss these two cases. 

As the first step, we shall write the general a(3^ integral given in Eq. (l2ip into the ILIs. 
With the standard manipulations, such as combining factors in the denominator with the 
UVDP parametrization and making translation of loop momenta, we can then get rid of the 
cross terms of momenta in the denominator. Some useful formula and further discussions on 
the UVDP parametrization method are given in Appendix B, which enables us to reexpress 
Eq. (l2Tl) into the following form 
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1 1 1 

(l+'Wl)"- 1 (l+v 2 y- 1 (1+V3P- 1 



{^(k 2 - m 2 ) + ^(k 2 - m 2 ) + j^K-h -k 2 + pf - ml]Y+^ 
T{a + + 7) f°° -A dvi A 1 f d 4 k[ f d 4 k' 2 
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,2 J 



Lv 1+^1 ^ 1+1)3 ' 1 (2+di+-u 2 )(1+u 3 )' v 2 ^ 3 + „ 1+t , 2+ „ 3 ^ Z^j=l l +Uj .J 

(22) 

where we have used (i=l,2,3) to denote a, /3, 7 respectively, and made the following 
momentum translation 

k'x = fci + ^ + ^ (p-k 2 ), (23) 

2 + f 1 + v 3 

k'2 = h+ 1 + Va P- (24) 

3 + f 1 + v 2 + v 3 

Below we shall drop the prime on k'i for simplicity. It is seen that the cross term of 
momentum is eliminated. 

In general, the power indices a, (3, 7 are positive integers, so that we have a + (3 + 7 > 3. 
Thus the above integral is convergent with respect to one of loop momentum fcjS. From the 
general structure of Eq. (|2ip . it is clear that the final result is independent of the integration 
order over k\ and k 2 . Without lost of generality, we can first integrate over k\ and explicitly 
obtain the expression: 
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(25) 



with 



m2 = (2 + v 1 + v 3 )(l + v 2 ) (2 + v 1 + v 3 )(l + v 2 ) A m) 

(3 + v x + v 2 + v 3 ) 2 P 3 + v x + v 2 + v 3 f^l + vj' 
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p( , (2 + Vi + v 3 ) a+ ^-\l + ViY^jl + ^ 3 ) 1 " 7 (1 + ^) a+7 ~ 3 m 
{Vk) - (3 + v l + v 2 +v 3 )^- 2 ■ {Zf) 



The above result is symmetric under the interchange between v\ (mi,a,i = a) and 
V3 (777.3, a 3 = 7)- 111 fact, the original expression is also symmetric under the permutations 
among v\ (m^aj = a), v 2 {m 2l a 2 = 0) and 7j 3 (7773,a3 = 7). Making the following scaling 
transformation for the momentum 



kl = (2 + ^ + ^ 3)(1 + ^^, (28) 
2 3 + v x + v 2 + w 3 



we then obtain the following more symmetric expression 
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with 



(27r) 4 [/2 - .M 2 (p 2 , m 2 ? ^)]a+/3 +7 -2 > 
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F{v h ) 



(3 + f 1 + t; 2 + v 3 ) 



2 



M * = r(-?L)-- — x - ?\ 

j^yi + vj 3 + Vl + v 2 + tj 3 

In the subsequent sections, we will show how this formula can naturally be obtained when 
merging the UVDP parametrization and the evaluation of ILIs with the Bjorken-Drell's 
circuit analogy. 

To go further, we need to consider some explicit values of a, ft, 7. As mentioned above, the 
only cases involving overlapping divergences are (1) a+/3+7 = 4 and (2) a+/3+7 = 3. Up to 
the field redefinition, we can always take the corresponding cases to be (1) a = 7 = 1, (3 = 2 
and (2) a = j5 = 7 = 1. We will consider these two cases separately in detail. 



IV. TREATMENT OF DIVERGENCES IN THE UVDP PARAMETER SPACE 

Let us first consider the simpler case with a = 7 = 1, = 2, where the general form of 
a/?7 integral Eq. (|25p can be simplified into the following form 



! 121 
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I n ^ - e tt^ } (3 + Vl + V2 + t, 3 ) 2 (i + v 2 ) 
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{2Tty[kl-M{p^mlv k )Y 
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1 f°° v— > 1 

;i6^ / II dv ^ 1 - E ! + ' (3 + «j + t; 2 + t; 3 ) 2 (1 + v 2 ) 

% — 1 J — 1 

M 2 M 2 

(30) 

The integral over the loop momentum k 2 above is logarithmically divergent, which represents 
the overall divergence. When carrying out the integration over the loop momentum k 2 , we 
have applied the LORE method to regularize it, where M. 2 = Ai 2 + /x 2 with /i 2 playing 



the role of IR divergence cut-off. In the following, we are always working in the massive 
theory, so there is no IR problem and the scale [x s can safely be set to [i 2 = 0. Also, in the 
limit M 2 oo, yi (&) ~ 0(^f ) 0, so yiS vanish identically. By power counting, the 
only contribution from y^s may arise when the overall quadratic divergence is multiplied by 
2/o( jtfi)- Such a contribution is finite and does not disturb the divergent terms. Nevertheless, 
in the present paper, we only focus on the divergent part to show the consistency of the 
LORE method. Thus, the limit M 2 — > oo is always taken and all the terms yi(jjz) are 
dropped below to simplify our expressions. 

It is not difficult to see that there exists a divergence in the region of UVDP parameter 
space with V\,v% — >■ oo, which reflects the divergence of subdiagram cry. To extract the 
divergence, we may focus on the region where V\, V3 > V with V ^> 1. In such a region, v 2 
0, which is ensured by the delta function, so the domain of the integration is transformed 
into dvi §y dv% J °° dv 2 . With such a treatment, it is easy to check that M. — > m\ and 
some terms small in comparison with V\ and v$ can be neglected. Thus the integral /121 is 
simplified into the following form 



'121 
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1 M 2 f°° 1 

= -7T7T^ ln ^-^) / rf ^^17' ( 31 ) 

(167T 2 ) 2 Wlj Jy V 1 + V 

where we integrate over v 2 and V3 as they are convergent. The remaining integration over v\ 
is divergent, which has to be regularized appropriately. The LORE method has been shown 
to be more suitable in this situation jl, 0], because such a divergence is a kind of scalar type 
divergent ILIs, which is the object that can be regularized in the LORE method, rather than 
other physical objects, such as propagators or the dimension of the theory. To regularize 
the UVDP parameter integral, it is more useful to transform it into a manifest ILL For that, 
one just needs to multiply a free mass-squared scale q 2 to vi, which will be determined by 
a suitable criterion. Eventually, it will cause the harmful divergences of different diagrams 
to be canceled. In general, such a scale can be the function of the intrinsic quantities in the 
theory, such as masses of particles or external momenta. 

1 f°° 1 



q 2 vi + q 2 V J {q 2 )v q{ + q 2 V 



M 2 

= ln^-7„, (32) 

where we have defiened q\ = q 2 v\. In the following, we will frequently encounter similar 
divergent integrals in the UVDP parameter space. Unless specified explicitly, we shall always 
use this prescription to deal with them. Here we would like to emphasize that the above 
prescription is the only one consistent within the framework of the LORE method. For 
other regularization schemes, such an approach cannot give consistent results, such as the 
Pauli-Villars regularization in which the regularized objects are the propagators of internal 
particles. 

It can be shown that in other regions of parameter space, there are no further divergences. 
Namely, besides the overlapping divergence given above, they contains only the harmless 



overall divergence from the integration of hi- Thus the general form of overlapping divergence 
in the integral 1\2\ can be written as: 



1 M 2 M 2 

In order to show the exact cancelation of harmful divergences for I±2i, it is necessary to 
calculate its corresponding counterterm diagram (cry), (see Fig. ([Zj)) 

r (c)(a 7 ) _ f d A k 2 1 , f d 4 k t 1 
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77? 



where DP{} denotes the divergent part. Such a counterterm integral can be easily computed, 

41 =+7T^( ln ^|-^)-(ln^-7.) (35) 

where the second factor comes from the subintegral (cry) part contained in DP{} and the 
first one from the integration of internal loop momentum hi- 

It is obvious that there is an exact correspondence between the factors in each expression. 
When taking the free scale to be /i 2 = 2q 2 V, the two divergent terms cancel each other 
exactly. Here the divergence contained in the UVDP parameter space in the region Vi,v 3 — > 
oo reproduces that of subintegral (cry), namely the integration over k\. We also notice that 
the divergences of J121 are factorizable and can be written as the product of two divergent 
integrals, i.e., one from the integral for the overall divergence and the other from the 
subintegral h\ (07) for the sub- divergence, which is transformed into and represented in 
the UVDP parameter integral of the region fi,f3 — > 00. This is the general feature when 
using the LORE method to disentangle the two-loop overlapping divergences. We will make 
a more explicit demonstration on this feature below by merging with the Bjorken-Drell's 
analogy between Feynman diagrams and electrical circuits. 



V. EVALUATION OF ILIS AND BJORKEN-DRELL'S ANALOGY BETWEEN 
FEYNMAN DIAGRAMS AND ELECTRICAL CIRCUIT DIAGRAMS 



In order to generalize the correspondence between the divergences in the UVDP parameter 
space and those in the subintegrals to more complicated cases, it is interesting to observe 
that the UVDP parametrization and the evaluation of ILIs in the LORE method naturally 
merge with the Bjorken-Drell's analogy between the general Feynman diagrams and the 
electrical circuits. A detailed description for such the analogy is referred to the book by 
Bjorken and Drell[l5j. It was originally motivated for discussing the analyticity properties 
of Feynman diagrams from the causality requirement. Here let us first establish such the 
analogy by developing a standard procedure and notation following Bjorken and Drell, and 
then apply it to the general a/?7 integrals by merging it to the LORE method. 

For a general connected Feynman diagram, we shall always denote the external momenta 
of the diagram by pi, ...,p m with the direction of coming into the diagram. Thus, according 
to overall momentum conservation, we have: 



in 
8=1 



(36) 



To each internal line we assign a momentum kj with a specified direction and a mass rrij. 
At each vertex, we have a law of momentum conservation of the form 



3=1 8=1 

where is chosen to be +1 if internal line j enters vertex i, while —1 if internal line j leaves 
vertex i, otherwise is defined to be 0. <fj S has the similar definition for the external lines 
which, by convention, are always taken to enter vertices. 

Each diagram has a definite number k of internal loops. However, we have the freedom 
to choose the concrete internal loops and assign each loop a momentum l r which are going 
to be integrated out along the loop. Thus, for each internal line j, we have the following 
decomposition: 

k 

kj = qj + 2J Vjrlr, (38) 

r=l 

where r\j r is chosen to be 1 if the jth internal line lies on the rth loop and the momenta kj 
and l r are parallel, and -1 if the jth line lies on the rth loop but kj and l r are antiparallel, 
otherwise r]j r is 0. Notice that here we introduce another kind of internal momentum qj, 
which will be determined after we adopt the UVDP parametrization for combining denom- 
inators to evaluate the ILIs. From the decomposition Eq.f l38|) . we can immediately obtain 
the following momentum conservation law for each vertex in terms of qf 

n m 

^ Ujqj + UsPs = 0, (39) 

3=1 8=1 

which follows from Eq.( l3Tj) and 

n 

= 0, (40) 

0=1 

which is a consequence of the definitions of and r]j r given in Eqs. ( 1371) and ( 1381) . 
The general structure of the Feynman integral can be written as follows: 

f N 

:,Pm) = J dh...d h ^^^ ^^^^ (41) 

where N represents the numerator of a general matrix element, which can be the products of 
external momenta, internal momenta, spin matrices, wave functions and so on. By adopting 
the UVDP parametrization, the above integral can be written as: 

u \ [*i ^ r (g=i%) rfr dvj x „ A 1 , 

I(p 1 ,...,p m ) = / d%...d% J II ttt w^i ^ 1 ~ T^~) 

J r(ai)...r(a n ) J A 1 (1 + „.)a.+i 1 + v . 



N 



LZ^3=i \ +Vj J 



/A..Afn^d-tr 



T( ai )...T(a n ) J i '" Jq A1 (l + ^) ai+1 v ^l + v 

N 



i=i y ' '"' 3=1 ' ~ J 



In order to obtain the required ILIs, we need to eliminate the cross terms in the denominator 
which implies 

V^^ = (43) 

3=1 J 

for each loop r = 1, ...,k. Now we have the enough conditions Eqs. fl39l) and f|4*3l) to 
determine the momenta qj for each diagram. The above procedure is essentially equivalent 
to the usual way of shifting the loop momenta for completing the square in the denominator. 
Our next task is to diagonalize the momentum integration variables so that we can integrate 
over each momentum integrals separately in Eq. ( |42lh 

Before doing the calculation, let us try to understand Eqs. ( |39|) and f T4"3l from an alterna- 
tive interesting perspective. First we put them into a more heuristic form: 

T T^- = 0, (44) 
• , l + v j 

qj in loop r J 

(Qj+Ps) = 0, (45) 

qj, p a entering vertex i 

we then arrive at a complete analogy between the Feynman diagrams and electrical circuits. 






FIG. 8: Left: Current conservation for each vertex i; Right: Conservation of the voltage in any 
loop r 



Specifically, we can think of the Feynman diagram as an electrical circuit and associate the 
momenta with the currents. Thus qj are the internal currents flowing along the circuit and p s 
the external currents entering it. When associating the parameters — ^y- with the resistance 

of the jth line (so Vj can be regarded as the conductance of the jth line), we see explicitly 
that Eqs. flH|) and (j4"5]) simply become the Kirchhoff 's laws in this circuit analogy. Eq. fj44"|) 
shows that the sum of "voltage drop" around any closed loop is zero, and Eq. fj4"5l) indicates 
that the sum of "currents" flowing a vertex is zero. 

Moreover, if we associate the voltage with the coordinate x M of the vertex, we can even 
inquire the physical meaning of Ohm's law: 

V = IR (46) 

to be the following relation by translating it into the language of Feynman diagrams: 



where qj, Vj are the momentum and UVDP parameter carried on the internal line, and 
Axj is the displacement between two points connected by this line. In fact, Eq. ( 14 7p is just 




V 



R 



FIG. 9: Ohm's Law 

the equation of motion for a free particle, which becomes more apparent in terms of the 
component forms, 



At 4 



1 



1 • i'j ■■ l + Vj l\tj q 3 



At,- 



1 

o ' 



(48) 



As the parameter is positive definite, the causal propagation of the particle is guaranteed: 

^w- = > 0- (49) 

As the particle goes in the qj direction according to ( 148]) , it moves either forward or backward 
in time depending on whether the sign of the energy qf is positive or negative. This agrees 
with the interpretation of causality of Feynman propagator in QFT. 

The above description provides us a physical picture of the circuit analogy which can be 
summarized as follows 



Feynman diagrams 


4$ 


Electrical Circuit diagrams 


(50) 


Displacement Axj 




Voltage 


(51) 


UVDP parameter 




Conductance > 


(52) 


Free particle equation of motion 




Ohm's law 


(53) 


Cross term cancelation condition for ILIs 




Kirchhof f's Law 


(54) 



while the positivity of the UVDP parameter Vi as the "conductance" is related to the causal- 
ity of propagation for the free particles. 




FIG. 10: Analogue of Feynman diagrams and electrical circuit diagrams 



In order to carry out the integral over l r in Eq. (|42p . it is useful to make the quadratic 
terms of the momentum l r diagonal. First write it in terms of the matrix form 
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j,r,r' 



1 + Vj 



^/ r M rr /Z r / = L T ML, 
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M, 



EVjrVjr' 
l + Vj 

1 J 



(55) 



where L T = (/i, • • • ,l k ) is the transpose of the vector L and M„i is a symmetric matrix. We 
then diagonalize the matrix M by an orthogonal transformation O with 

L = OL', T MO = diag(Xx,- ■ ■ ,X k ) = diag(\ + , ■ ■ ■ ,A_( fc _i)), (56) 

where A r (r = 1, • • • , k) or A + , A_( r ) (r = 1, • • • , — 1) are the eigenvalues of the ma- 
trix M, corresponding to the eigenvectors V — ■ ■ • , l' k ) T = l'_n)i ' ' ' > ^'_(fc_i)) T - As 
the transformation matrix O is orthogonal, the integration measure remains unchanged 
dH[ ■ ■ ■ dH' k = d 4 h ■ ■ ■ dH k . Thus, the integral Eq.(fl2l) can be simplified to: 

/fe,..., Pm ) = r(ai) rK) y o ntTT^^-Ei^) 
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£(g^i) rfV A 1 

r( Ql )...r(a n ) 7 f = {(i + ^)- +l1 ^i + V 

/ rf 4 /;...^ s 2 (57) 



For a generic k-loop integrals where k > 2 and n > k, we have the inequality 57J™ =1 «j > 
fc( - fc 2 t " 1 ' 1 >2k — I. 1 Thus we can explicitly integrate out the loop momenta, as these integrals 
are already convergent. In particular, when the numerator N contains no l[ terms, we can 
integrate out the last {k — 1) internal loop momenta, say l' 2 ,l' 3 , ■ ■■,l' k - 



I(Pl,-,Pm) = 



k-1 



d 4 Z' = — 5 . (58) 

By a rescaling Z + — > y^X^l' + , we then obtain the following form: 
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with the definition of the determinant for the matrix M 

fe fc-i 



det |M| = JjA r = A+ JjA_ (r) . (60) 



r=l r=l 



The first inequality comes from the fact that in order that the expression is generic, we need to con- 
sider every type of internal momentum combinations in the denominator, such as — Pj) 2 , {l\ + — 
Pi(i+i)) 2 , (^ + ^2 +—+^fc — Pi2...k) 2 ■ The total number of the combinations is k+{k— l) + ... + l = M*j+ii. 
If all the types of combinations appear in the denominator, then the inequality holds. 



The above expression is the required form of ILIs, where the ILIs for the momentum integral 
on l + reflects the overall divergence of the Feynman diagram. From the above expression, 
it is clear that the UV divergences contained in the loop momentum integrals over l'_r r ) 
(r — 1, • • • ,k — 1) for the original loop subdiagrams are now characterized by the possible 
zero eigenvalues A_( r ) — > (r = 1, • • • , k — 1) of the matrix M in the allowed regions of 
the parameters Vi (i = 1, • • • , n). Namely, each zero eigenvalue A_( r ) — > resulted from 
some infinities of parameters Vi in the UVDP parameter space leads to a singularity for the 
parameter integrals, which corresponds to the divergence of subdiagram in the relevant loop 
momentum integral. 

By applying the general LORE formulae to the above integration over the momentum 
l + , we get: 

T( \ = r (g = i^-2A: + 2) r°°" dVi " 1 1 



N 

lim V / d% 



i=l ^ " j=l 

i(_l)E"=i«i 



N,M? ' / [V^n gj~ m 

r(E- = i«, -2^ + 2) r°°" dVt r/l ^ i 



r(« 1 )...r( 



2k + 2 ) rfr d Vi ^ i 



(det |M|) 5 
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with 



a = ^«i-2fc, ^ 2 = ^K 2 -g|)/(l + t;,) (62) 
j=i j=i 

where I^ 2a (.-M 2 ) is the regularized 1-fold ILI for the possible overall divergence of the Feyn- 
man diagram. 

In general, there are (k — 1) zero eigenvalues A_( r ) — > (r = 1, • • • , k — 1) in the UVDP 
parameter space for the k-rank matrix M, and they correspond to the divergences of the 
(k — 1) loop subdiagrams in the momentum space. In principle, to arrive at the k-fold ILIs 
for the k-loop Feynman diagrams, one may perform (n — k — 1) integrations in the UVDP 
parameter space. It requires us to appropriately analyze the zero eigenvalues of the matrix M 
associated with the corresponding regions of the UVDP parameters. Alternatively, one may 
make an appropriate parameter transformation, so that the integrations on the (n — k — I) 
parameters become convergent for the considered regions of parameters in a new UVDP 
parameter space, thus they can be integrated safely. As a consequence, we obtain the 
desired k-fold ILIs. We shall illustrate in detail its consistency and advantage by applying 
it to the general af5 r j integrals in the 4 scalar theory. So far it becomes apparent that the 
above general procedure explicitly realizes the UVDP parametrization and systematically 
obtains the ILIs and this shows the powerful advantage when merging the LORE method 
with the Bjorken-Drell analogy between Feynman diagrams and electrical circuit diagrams. 

In order to demonstrate explicitly the correspondence between two kinds of divergences 
in the UVDP parameter space and in the momentum space, we are going to apply the above 
general procedure to the a/?7 integral in next section. 



VI. DIVERGENCE CORRESPONDENCE BETWEEN SUBDIAGRAMS AND 
UVDP PARAMETERS 



The corresponding Feynman diagram for aftj integral is already shown in Fig. (jSJ). With 
the internal momenta kj and the particular choice of loops defined therein, we can rewrite 
the a/?7 integral as follows 
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where we have introduced a new notation ctj (i=l,2,3) corresponding to a, /3, 7 in the second 
line. According to the diagram, we have the momentum conservation, either for overall 
diagram or for both vertices: 

Pi = -P2 = Pi 

and 

p l — k 1 — k 2 — k 3 — 0, (64a) 
p 2 + k 1 + k 2 + k 3 = 0. (64b) 

Following Eq. fl38|) . we decompose the internal momenta kj into two parts: one represents 
the loop momentum flowing along the line j, and the other for the external one carried by 
3 

h = qi+h, (65a) 
k 2 = q2 + h-, (65b) 
h = qs-h-k, (65c) 

We then arrive at the momentum conservation laws for either vertex in terms of qj 

Pi = qi + <?2 + <?3 = -P2 = P- (66) 
Replacing the kj with qj and l r in Eq. (|63|) and changing the integral variables to l r give us 

f d% r d% r(a + /3 + 7 ) rA dv * m V 1 
a ^ J (2vr) 4 J (2tt) 4 r(a)rG8)r( 7 ) J 11 (1 + ~ i + 



with 

3 „2 _2 

V\ 1 + ^3 1+^2 1 + ^3 



D = Y q -^-^ + 2(-^ ^-)J a + 2(-^ ^)/ 2 + L T ML, (67) 



where we have introduced the definitions: 



The elimination of the cross terms in the denominator D requires that 



T— T - r = 0, (68a) 

1 + V i 1 + f 3 

' /2 0. (68b) 



1 + V 2 1 +V 3 

These two formula explicitly illustrate the Kirchhoff's law for two loops in the electrical 
circuit analogy of . By taking into account Eqs. ( 1681) and ( )66l) together, we obtain the 
solutions: 

1 + Vl (an \ 

Qi = : : P, 69a 

3 + Vi + V 2 + V 3 

<?2 = =- : — V, 69b 

3 + v 1 + v 2 + v 3 

93 = =- : : p. (69c) 

3 + V 1 + f 2 + v 3 

In order to perform the integral over l r , we may first diagonalize the matrix M by a 2 x 2 
orthogonal matrix transformation O, so that 

L = OL', T MO = ( X ' x ° 

with Ai 2 = A +> _ given by 



(l + T +-)± A /(l + T +-) 2 + 4A 

A ± = W V 2 1+ " (70) 

A = d6t|M| = (1 + vMl + v 2 ) + (1 + v 2 )(l + v 3 ) + (1 + v 3 )(l + v,) ' (71) 

which are the two eigenvalues of the matrix M corresponding to two eigenvectors V = 
V 2 ). Since the transformation matrix O is orthogonal, the integration measure remains 
the same dH^dH^ = dH\dH 2 . Thus, the a/3 7 integral can be reexpressed as: 

3 



, d% r d% rja + p+j) f°°yr dvi v x 



J {2-ny r(a)rC9)r( 7 ) 7 l\ (1 + Vi )^ v f^i + 
i 



'4/' ^ ^/4;/ tV„ i /O i _,\ roo 3 



d% f dH'_ I> + /3 + 7 ) 



(2tt)*./ (2vr)4r(a)r(/3)r( 7 )7 l±(l + «,)«*+! v tf 1 + ^' 
1 



73+7 



(72) 



From this formalism, it can be shown that the integration over l'_ represents the integral 
over subdiagrams, while the one over l' + is an overall diagram. However, the matrix M is 
not always invertible, since the determinant of M vanishes when any two of ViS tend to 00. 



More specifically, take t>i,t>3 — > oo for example. In this case the eigenvalue A_ vanishes. 
It is also noted that the combination and A_ZL in Eq. (1721) are on equal footing in 

the denominator, it is then expected that and A_/'_ approach to infinity at the same 
speed when both l' ± — > oo. Thus, when considering A_ — > while keeping A + finite, it 
requires that the speed of l'_ tending to infinity is faster than that of l' + in order to keep 
the balance. Recall that in our previous general discussion on the divergence behavior of 
overlapping diagrams, one of the features for the sub divergences is that the integration 
variables approach to infinity faster than the overall one. 

Based on the above analysis, we may conclude that the integral over l'_ reflects the 
asymptotic behavior of subintegrals when the corresponding UVDP parameters approach to 
infinity. Here we would like to emphasize that the integration over l'_ does not correspond to 
any particular loop in the original Feynman diagram. Rather, it represents all subintegrals 
and is specified according to the asymptotic regions in the UVDP parameter space. For 
instance, when the divergences in the UVDP parameter space occur in other regions, such 
as Vi, t>2 — > oo, then l'_ reflects the loop composing of lines 1 and 2. The above explicit con- 
struction helps us to understand the intuitive analogy between the Feynman diagrams and 
electrical circuits. Especially, it illustrates why and how the divergences in the subdiagrams 
are transmitted to the corresponding divergences in the UVDP parameter space. Let us 
further demonstrate this point from another perspective. By explicitly integrating over l'_, 
we obtain, 



'a/37 



d% r(q + P + 7 ~ 2) 

16tt 2 J (2tt) 4 r(«)r(/3)r(7) 

1 

Ei=i 45? + A+/?r^- 2 



1=1 v ' J=l J 



(73) 



which explicitly shows that when A_ goes to zero, that is, any two of the three UVDP 
parameters V{ approach to infinity, the integrand becomes singular and the integrations over 
the UVDP parameters give some UV divergences. 
By defining a new integral loop momenta Z + as 



l + = y/X + l' + , (74) 
we can transform the a/?7 integral into a more tractable form: 



'a/37 



i f d% r(a + /3 + 7 -2) f°° yr dvj V 1 ' 

i6tt 2 7 (2vr) 4 r(a)rO0)r( 7 ) J f = {(i + ^) Qi+1 + 



1 



(det |M|) 2 [V 3 Q]-m] /2 lQ+/3+7 _ 2 
i r(a + /3 + 7-2) PA d Vi A _}_.„, , 

^ r(a)r(/3)r( 7 ) J f = J- (i + 1 ^i + V w 



16vr 
dH 



(2tt) 4 [11 - M\ V \ ml, v k )}«+^~ 2 



(75) 



where 

(i + Vl ) i - a {i + v 2 ) i ~p{i + v 3 y-' 



F{Vk) ~ (det\M\Y 

(l + Vl ) 3 ~ a (l + v 2 ) 3 ~^l + v 3 ) 3 ^ 



(3 + vi + v 2 + v 3 y 

™? 

+ Vj 3 + v\ + v 2 + 



3 2 2 3 2 



i=i J i=i 

The above equation is equivalent to the form in Eq. (!25|) with a rescaling given in Eq.( |28|) . 
Nevertheless, the derivation here is more general and systematic and it explicitly shows the 
advantage when merging the UVDP parametrization and the evaluation of ILIs with the 
Bjorken-Drell's electrical circuit analogy of the Feynman diagrams. 

In general, the integration for the momentum l + is divergent and needs to be regularized. 
By applying the LORE method to the momentum integral, we obtain: 



'a/37 



I> + /3 + 7 -2) f -A dvi 



Jo f =1 ( 1 + v ^ ~t 1 + v i 



i6tt 2 r(a)r(/5)r(7) 
- r d 4 i 

lim Vr N / — — - 

N,Mff^ 1 J (27T) 4 [/2+M2+7W2( 2 2 ) ]a+/3+7 -2 



r ^ N f d% i{-iy+^ 

L^l 1 / (0^\A \J2 — I — M 2 4- A/f2fr>2 th 2 „, ^la+/3+7-2 V ' 



-^ r( ° + .f.; 7 r - 2) r n - £ ^ w 

i6vr 2 r(a)r(/3)r(7) j f = |(i + ^) 2 ^i + ^i 



(a+/3+7-4) 



r- 

When applying the above general formula to the case a = 7 = 1, = 2, with the similar 
calculation as the one in the previous section, the result is the same due to the equivalence 
of Eq. Q25) and Eq. (j75]h 



'121 



A dvj A 1 (l+^) 2 (l + t; 2 )(l + t;3) 2 * 

I — 1 J — 1 



-1 poo 3 3 

l_1 nw-xd 



1 



167T 2 ) 2 Jo A _! 1 + ^' (3 + + ^2 + ^ 3 ) 2 (i + U 2 ) 

M 2 M 2 

+*Ajjt» < 78 > 

which shows that the singular behavior in the region 1)1,1)3 — » 00 becomes obvious as 
det |M| = A = due to the zero eigenvalue A_ — > 0. For the other two regions: v±, v 2 — > 00 
and v 2 ,V3 — > 00, the additional factor Tjj^ in these two cases makes the integration finite. 
In contrast, for the case a = /3 = 7 = 1, there is no such a factor, so that there are more 
UV divergent structures in all the three regions andis going to be discussed in detail below. 



VII. TREATMENT OF OVERLAPPING DIVERGENCE AND ADVANTAGE OF 
THE LORE METHOD MERGING WITH BJORKEN-DRELL'S ANALOGY 

This section shows that the LORE method merging with the Bjorken-Drell's circuit anal- 
ogy has the advantage in analyzing the more complicated and challenging overlapping di- 
vergence structures of Feynman diagrams. For an explicit demonstration, we are going to 



consider the case with a = /3 = 7 = lin the a/?7 integral. The difficulty lies not only in 
the quadratic divergence but also in the more complicated overlapping divergence structure. 
It will be seen that the LORE method merging with Bjorken-Drell's analogy is extremely 
powerful in unraveling the overlapping divergences. 

The general form of a/?7 integral (Eq fTSl) can be simplified to: 



% f°° -A- dvi 1 1 

111 = 16^J 11(1 + ^)2^ _ j^T+^'Jd^. 

N 

lim y a 

M A/r2 Z — / 



\M\y 



v - dH, 



N,Mf 
' 1 1=1 



W ELi tzut + P+ + M? 



(16t 2 ) 2 / n^ + Vi) 251 - 1 5 1 + V(3 + »i+»'2 + «3) 2 

I 1 J 1 

Wl-M\\n^ 2 - luJ + l)l (79) 

where we have regularized the overall quadratic divergence of loop momentum integral by the 
LORE method. The mass factor Ai is given in Eq.(??). The UVDP parameter integrals are 
more involved due to the appearance of the overlapping divergences. From the expression of 
integral Jm, it is seen that the three subintegrals 07, (3j, and a(3 are all divergent. With the 
analogy of circuits, we have shown that the UV divergences arising from the large internal 
loop momenta transmit to the asymptotic regions of UVDP parameter space, where the 
divergent conductances correspond to the following asymptotic regions in the circuits: 

subdivergence in aftj diagrams divergence in UVDP parameter space 

Circuit 1 : cry divergence v± — > 00, v% — > 00, v 2 — > 0, (80) 

Circuit 2 : f3 r y divergence <^ v 2 — > 00, v 3 — > 00, V\ — » 0, (81) 

Circuit 3 : a(3 divergence <=>■ V\ — » 00, t> 2 — > 00, V3 — > 0. (82) 

This result can also be obtained by considering the singularities in the determinant det \M\ = 
A as discussed in the previous section. 

Note that Eq.f J79|) has a permutation Z3 symmetry among the three pairs of parameters 
(vi,mi), (v2,m 2 ), (t>3,m 3 ), so the treatment on three asymptotic regions in the circuits is 
essentially the same. Let us consider in detail the first case of the Circuit 1. 

Circuit 1: 17 — > 00, V3 — > 00 and v 2 — > 0. In such region, the integral domain can be 
written as jy dv± jy dv% with M 2 — > m\ and F(vj) — > ^ 1+ ^_^^i 3 ^ ■ Thus the integration 
is simplified to: 

r ( )(a 7 ) _ 1 r dv ^ r dv * (i+^i) 2 (i+^ 3 ) 
-11 



4 1 , ir; ,n • o s2 



Ml 



[M c 2 -m^(ln^--7a, + l)] 



m 



2 



00 ^ 

dv\ — + .. 



(16tt 2 ) 2L c zy ml n J v v x + V 

_L_|M 2 - „ 2 (l n § 2 - lu + l)](ln^L - 7 „) + .... (83) 



Note that in the last step, we have applied the LORE method with the treatment discussed 
in Eq. (l32|) . The dots represent other terms, such as single logarithmic divergent term and 
finite terms, which are irrelevant to our discussions as our main purpose here is to check the 
cancelation of the harmful divergences. Note that our result here is factorizable. 

In order to compare the above divergence structure with those contained in the subdia- 
gram (cry), we calculate the counterterm diagram I^q : 

r (c)(a7) _ / ^ fc 2 1 Dpr f d%_ 1 ; 

111 J (2tt) 4 k\ - m 2 2 V (27r) 4 (Jfe?-mf)[Jfe§-m§]' 

1 M 2 M 2 

-(i^T (ln /" 7 - )[M ^^ (ln 4- 7 - +1)] ' (84) 

where DP{...} means the divergence part of the integral in the bracket, and fi 2 is the 
renormalization scale. It is then manifest that when we choose /x 2 = 2q 2 V, the harmful 
divergence parts in the two expressions cancel exactly. 

With a similar discussion based on the permutation Z3 symmetry, it is easy to show that 
the harmful divergent parts in the Circuit 2 and Circuit 3 also cancel exactly. So far we 
prove that there is no harmful divergence for the case a = /3 = 7 = 1 when combining with 
the corresponding counterterm diagrams. We would like to mention that this is different 
from the dimensional regularization, that in that we have an extra term corresponding to the 
quadratic divergence M 2 . As emphasized in [3j, this term is natural to maintain the correct 
divergent behavior of the original diagram, which can play an important role in effective 
field theory for obtaining the correct gap equation to describe the dynamically generated 
spontaneous symmetry breaking (Io|. It will explicitly be shown below that the presence of 
this term prevents us from having a mass independent renormalization scheme. Thus, a 
consistent renormalization with a well-defined subtraction scheme must be proposed for the 
LORE method. We propose the following subtraction scheme: 

(i) For quadratic divergence [M 2 — M 2 ), subtract [M 2 — fi 2 ) and leave (/1 2 — M 2 ) in the 
finite expression; 

(ii) For logarithmic divergence (log jfe — 7 W ), subtract (log —j- — 7u>) and leave term log j-^ 
in the finite expression. 

Such a scheme may be regarded as a kind of energy scale subtraction scheme at fi 2 and 
is similar to the usual momentum subtraction. For the logarithmic divergence, it appears 
to be a MS'Tike scheme in the dimensional regularization as it is associated with the Euler 
number 7^ = 7^. It is interesting to note that once the energy scale subtraction scheme for 
both the quadratic and logarithmic terms is set up at the one-loop level with a correlated 
form (M 2 — fi 2 ) and In M 2 //i 2 via a single subtracted energy scale fi 2 , and suppose that the 
correlated form with a single subtracted energy scale /1 2 is required to be maintained, thus 
either the rescaling fj 2 — > e a ° fi 2 or shifting /i 2 — > fi 2 — a§m 2 for the subtracted energy scale 
fi 2 will not be allowed. As a consequence, the mass renormalization at higher loop becomes 
well-defined through such an energy scale subtraction scheme at the one loop level, namely 
fixing the correlated form for the quadratic and logarithmic terms via a single subtracted 
energy scale fi 2 . 

Based on the above analysis and discussions, we arrive at the following theorems: 
Factorization Theorem for Overlapping Divergences: Overlapping divergences which con- 
tain divergences of subintegrals and overall one in the general Feynman loop integrals become 
factorizable in the corresponding asymptotic regions. 

Substraction Theorem for Overlapping Divergences: For general scalar-type two-loop in- 
tegral I a p~f, when combined with the corresponding subtraction integrals (which is composed 



of divergent subintegrals multiplied by an overall integral), the sum will only contain harm- 
less divergence. 

For completeness, we have the following theorems for dealing with the Feynman integrals 
which do not involve the overlapping divergence. They are so obvious that the proofs are 
omitted here. 

Harmless Divergence Theorem: If the general loop integral contains no divergent subinte- 
grals, then it is only possible to contain a harmless single divergence arising from the overall 
divergence. 

Trivial Convergence Theorem: If the general loop integral contains neither the overall 
divergence nor the divergent subintegrals, then it is convergent. 

In summary, the LORE method can properly deal with the overlapping divergences, es- 
pecially the sub divergences which is transformed appropriately into the divergences in the 
UVDP parameter space. To extract them, we need to explore the integrals in different 
asymptotic regions of the parameter space. Moreover, we demonstrate that these overlap- 
ping divergences can well be treated by the LORE method when merging with the Bjorken- 
Drell's analogy between general Feynman diagrams and electrical circuits, especially the 
correspondence between the UVDP parameters and the conductances of internal lines in 
the circuit analogy. By applying this intuitive picture, we can immediately recognize how a 
divergence in the region of UVDP parameter space corresponds to a certain original diver- 
gent subintegral composed by the lines that the divergent UVDP parameters are attached 
on. This correspondence also helps us to find the right counterterm diagram to cancel the 
notorious harmful divergences. As a result, we are left with only the finite terms and the 
harmless divergence which can be absorbed into the overall counterterm at two-loop order. 
These results are summarized in the four theorems presented above. It is interesting to 
note that the extension of the LORE method to the calculations beyond two-loop order is 
straightforward, although we are aware that the degree of complication and difficulty in- 
creases dramatically with the increase of loop orders as more and more Feynman diagrams 
are involved. Nevertheless, it is clearly indicated that merging with the Bjorken-Drell's anal- 
ogy between Feynman diagrams and electrical circuit diagrams, the LORE method gets its 
consistency and advantage in the multiloop calculations, especially with the aid of computer. 



VIII. APPLICATION TO ^ THEORY AT TWO-LOOP ORDER 

The discussion and analysis in the previous sections on the general two-loop integrals 
appear to be a little bit too abstract. In this section, we shall take the simple scalar 4 
theory as a concrete example to illustrate the LORE method in a practical calculation, and 
leave another application involving tensor-type integrals to a separate paper (22|. 

The Lagrangian density for </> 4 theory is: 

C = \d^d^ - \m 2 <? - (85) 



Its Feynman rules may be found in the standard textbooks, such as [16|, Il7j . Our main 
purpose here is to explicitly calculate the two-loop contributions to the mass term and 
coupling constant from the self-energy and vertex diagrams. From this practical calculation, 
we will demonstrate in detail the consistency and advantage of the LORE method when 
merging with the Bjorken-Drell's circuits analogy. 



A. Renormalization At One-Loop Level 



Before proceeding to a detailed calculation at two-loop level, we need the one-loop coun- 
terterms first. This is equivalent to specify the renormalization condition in the LORE 
method, which is the main goal of this subsection. 

At the one loop level, there are two types of diagrams corresponding to the self-energy 
correction and vertex correction respectively, as shown in Figs. (ITTj) and ( 1121) . 




P P 



FIG. 11: 




(s) (t) (u) 
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For the self-energy correction, the calculation is straightforward and the result is given by 

1 f d A k i 



(2tt) 4 k 2 - m 2 



^ r 2/, M r - / m 9X\n 

"W |Mc m(ln 4" 7 " + 1 + ! ' 2( 4 ))1 

-^ VM c- m2) - m2(la §-^ )] (86) 

where we have applied the LORE method to obtain the result in the second line. The 
result in the last line is obtained by taking /i s = and M c — > oo. Under the energy 
scale /i 2 subtraction scheme described in the previous section, the mass and wave function 
counterterms take the following forms: 

- *® = ^[(^ - - ™> f - 7«)], (87) 
ibf = 0, (88) 



and the finite term is found to be 



- iM '« = -2(S)5 [( ' ,2 - m2) - m2ln ^' (89) 



Note that the result given in Eq. fl89|) is different from the one obtained by using the dimen- 
sional regularization method. The difference arises from the quadratic behavior fi 2 in the 
renormalization counterterm. This difference may have important physical implications: it 



greatly changes the renormalization group |l2Ml4l| with a power law running, and generates 



the physically meaningful dynamical mass scales in the effective field theory 0, [Io| . 

By a similar calculation of the one-loop four-point Green function, we obtain the following 
vertex correction for s-channel: 



-iA 



(s) HA) 2 f d A k 



W " 2 J (2tt) 4 k 2 - m 2 (k + p) 2 - m 2 
X 2 f d A k f 1 , 1 



dx-, 



2 J (2vr) 4 J [xk 2 + (1 - x)(k + p) 2 - m 2 } 2 

X 2 f d A k f l , 1 
dx- 



2 J (2tt) 4 J [k 2 + x(l -x)p 2 -m 2 ] 2 

^S?L ^ m 2 - x{L x)p 2 ~ 7j ' (90) 

where — p 2 = — (pi + P2) 2 = s. For other two channels (t- and u-channels), we will obtain 
the same expression except for the definition of p 2 : —p 2 = —(pi —P4) 2 = t for t-channel and 
— p 2 = — (pi — P3) 2 = u for u-channel. 

According to the renormalization scheme of the LORE method, we have the following 
counterterm: 

.M) 3zA 2 ri M 2 

-^-w^" 1 " 1 (91) 

where the factor 3 comes from three diagrams corresponding to the s, t, u-channels. 
So the renormalized vertex correction is simply given by: 

iX 2 f 1 /i 2 /j, 2 /j, 2 

— iA(i\ — —. — — / dx\hi — „ ; ; — h hi — ; r — h In — ; —1.(92) 

(1) 2(4tt) 2 J 1 m 2 + x(l-x)s m 2 + x(l - x)t m 2 + x(l - x)u [ V ; 

B. Self-Energy Contribution at Two Loop 

There are two diagrams contributing to the two-loop self-energy corrections, which are 
shown in Figs. ([TBI and f JT4|) . 



The calculation of diagram (a) is straightforward and the result is: 

.,,2(a) 1/ . x s 2 f d * k l i i f rf4fc 2 



( 2 ) ' 4 V ; J (27r) 4 k\ -m 2 k\ - m 2 J (27r) 4 k 2 - m 2 
1 iX 2 M 2 . . „ 2/ M 2 

77T7T^ ln — \ ~ Tw) ' Wc ~ m 2 (^ - 7. + 1 • 93 
4 (Iott)" 1 m 2 m z 
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The corresponding counterterm diagrams are shown in (ai) and (02), and their sum gives: 
. , /f 2(ai)+(a2) 1 , . . s , s f d A k i 2 1 f , f d A k 



( 2 ) " 2 l n m V (2tt) 4 (P - m 2 ) 2 ^ 2 l a; J (2vr) 4 A; 2 - m 2 

-» l(i^ { ^ 2 - ^ - m2(ln f - ^ )](ln ^ ~ ^ 

M 2 M 2 
+ (ln -f - 7 J • [M c 2 - m 2 (ln -| - 7w + 1)]}, (94) 

where 5 m 2 is the one-loop mass counterterm defined in Eq. flH7|) and only the t-channel 
vertex counterterm, which is represented in Fig.l3(a 2 ). Thus, the sum of (a), (ai) and (a 2 ) 
gives: 

• , /f 2(a)+(al)+(o2) 1 ^ 2 r n \rfJ\/f2 2\ 2/i \i 

- lM (2) = I(I^j2" {-(lll 7l2--^)[( M c )~™ On -^2" -7*)] 



+[(// 2 -m 2 )-m 2 ln4]ln4}- (95) 



According to our renormalization scheme proposed in the previous subsection, the overall 
two-loop counterterm for diagram (a) is defined as: 

- = I ^ - 7.)[(M 2 - //) - m 2 (ln ^ - 7w) ], (96) 



and the renormalized correction to the two-point function from this diagram is: 

- = 777^ t(^ 2 - - 2 ) " - 2 ^ ^] In f^- (97) 



Let us now compute the most complicated diagrams in 4 theory at two loop order, 
namely the sunrise diagram. According to the internal momentum parameterizations shown 
in the diagram (b) and the general Feynman rules, we can write the expression explicitly: 



.,,2(6) 1/ -^2 f d 4 *7 f d 4 k 2 % % % 

lM (2\ = 7H A ) 



( 2 ) " 6 V ; 7 (2tt) 4 y (27r) 4 (A; 1 +p) 2 -m 2 A; 2 -m 2 (A; 1 + A; 2 ) 2 -m 2 



zA 2 /■ rf 4 A;! f d A k 
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6 7 (2tt) 4 y {2^ [{h + p) 2 - m 2 ]{k 2 - m 2 )^ + k 2 ) 2 - m 2 



(98) 
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Obviously, the above integral is just the special case of af3 r ) integral with a — f3 — j — 1 
and the same mass m 2 = m 2 . Thus we can apply the general result Eq. f!79|) to this case 



m _ /a 2 r A ^ _ A i n- =1 (i + ^) 2 
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with 



\ M c~M 2 (ln^-7. + l)], (99) 

.M 2 = m 2 p 2 . 

3 + V 1 + f 2 + v 3 

We now compute the contribution to the two-point Green function, rather than just 
giving the asymptotic expression for showing the cancelation of harmful divergences, which 
was already demonstrated in the previous section. For this purpose, it is useful to introduce 
a new set of UVDP parameters u,v,w via 



1 1 



l + <ui ' (l + u)(l+w)' 

1 _ u 

1 + V 2 1 + U ' 

1 1 



1 + W 3 ' (1+U)(1+V)' 

so the integration measure is now 
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(100b) 
(100c) 



''" '''' ;S(1-— -), (101) 



with 



(1 + ufJo (l + w) 2 (l + v) 2 K 1 + w l + v" 



M 2 = m 2 ; r P 2 , (102) 

(i + uY , x 

F = r , r 12 - (103) 

i M + (l+w)(l+v)\ 



With the above transformation, we finally arrive at: 

,,2(6) *A 2 f°° , [°° dw dv 1 1 



( 2 ) 6(16tt 2 ) 2 7 7 (1+ W )2(1 + ^) 2 ' 1 + " l + 

M 2 

[^-^ 2 (ln^-7« + l)]- (104) 
For the quadratic divergence, we can carry out the integration: 

_,m 2(6) - — m 2 r_*^ 1 1 * 
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6(16vr 2 ) 2 c 7 (1 + w) 2 (1 + w) 2 v 1 + w 1 + 
11 1 

o *''« + aT+T^ + I" + HT^imo] 2 (1 ~ (1 + »)(! + ») )] 

which is local when we choose the free scale g 2 = /i 2 . 

For the logarithmic divergence part, the result is given by: 

_.w2 (6 ) *a 2 r 00 rf w dt; i i v r , 1 + t* 

' (% 6(167r 2 ) 2 7 (1 + ^(1 + ^ l + 1 + M N+(TT^)] 2 

(1 + m)(1 + w){l + v)[u+ {1+w)(1+v) ] M 2 

From the analysis and discussion presented in previous sections for the overlapping diver- 
gences, there are three parameter regions which contain divergent contributions. To extract 
them, we need to separate the general expression into several parts, each of which may give 
an asymptotical result in a single region. In terms of the new set of UVDP parameters 
u,v,w, the situation becomes much simpler than the original parameters Vi,v 2 ,v 3 . The 
coefficients of the logarithmic divergence in the above expression can be separated into the 
following four parts: 

1 + U r 2 U 2 



- n ■ wn ■ : u ■ 1 



N + (i +w )W) ] (l + u)(l + w)(l + v)[u+ {1+w ] {1+v) ] 
-m 2 + - 5 — (1 — — -— z)m 2 



u + (1+4(1+^) ' i u + (i +w ](i + v) ] 2 " (1 + «')(! + '')■ 

V 



U 2 



= I + (11 + III) + IV. (107) 
These four parts are divergent in the following asymptotic UVDP parameter regions: 



I 
II 
III 
IV 



u — > oo vw = 1; v i — > oo v 3 — > oo v 2 — > 

f— t-OOM— 7-0 w— )■ 0; f 2 — )■ oo f 3 — > oo Vi — > 

w — > oo u — > v — > 0; ^2 — >■ oo v± — > oo v 3 — > 
p 2 > m 2 , 



where the second part in Eg. (jl07j) contains two asymptotic regions II and III symmetric 
under the exchange of parameters v and w or v% and v 3 . In general, it is difficult to carry 
out the whole integration and obtaining a complete result for —iMJ^ due to the complicity 

of Ai 2 in the logarithm. But it is sufficient for our present purpose to obtain —iM?^ 
by simplifying Ai 2 in the above asymptotic regions, allowing us to get the results up to 
logarithmic divergence. 

Region (I): u — > oo. In this region, Ai 2 can be simplified to the form: 

M 2 ~ m 2 . (108) 

So the approximate expression is given by: 

■ w2(&)(i) *A 2 f°° dw dv 1 1 s 



(2)W g( 167r 2) 2 y o (1 + M )2 (! + w )2 V 1 + w 1 + 



./o 



2 



u + 77- — w-, , , m 



m 2 (ln^|- 7w + l) 



'0 « T (l +tl) )(l + „) 

zA 2 m 2 M 2 M 2 
-^b(ln^-7. + 2)(ln^--7o, + l). (109) 

Notice that in the case u — > oo, the only relevant mass scale is m 2 , so we can simply take 
q 2 = m 2 . The final expression in this region is found to be: 



*«SgS? - -g(^" i2 i( ln § - + 3 ('» § - < 110 > 



Region (II+III): u — )■ oo or w — )■ oo. It is interesting to note that the integral is symmetric 
under the exchange of parameters v and w. Thus taking the limit v — » oo, or w — > oo, we 
shall arrive at the same results. In both cases, the asymptotic form of Ai 2 is given by: 

M 2 ~m 2 (111) 



2(b)(ii+ni) iX 2 [°° dw dv 



6(16tt 2 ) 2 J (l+w) 2 (l + v) 2 K 1 + w 1 + v' 
1 r du M[_ 

[1 (i + -)(i + -) J 7o [-+(tt4ttio] 2[ - 2 7w+ J 

i\ 2 2 Z" 00 dw dv r/1 1 1 , 



6(16tt 2 ) 2 7 (1 + w) 2 (1 + w) 2 v 1 + w 1 + v' 

[1 " 77- ^r- J (l + «0(1 + «)0n ^ - 7- + 1) 
(1 + + 

zA 2 2r f° rfw A 00 11/n M 2 



6(167r 2 ) 2 Jq 1 + w 7 1 + 1> m 2 

iA 2 2r „ M c x ,„ M r x , X 

"eOo^F" 1 [2(ln ^f -^)-l](ln-r|-7. + l), (H2) 

where in the third equality, we have used the constraint in the delta function (l + w)(l + v) = 
(1 + v) + (1 + w) to simplify the integral into the form of two 1-loop ILIs which can be 



regularized by the LORE method as shown previously. Again the only mass scale in the 
limit v — > oo or w — > oo is the mass of the particle m 2 , so the scale q 2 can be fixed to be m 2 
and the result is given by: 

- «4«" + "" K -^m>Wn § - + (In § - ■*.)]. (113) 

Note that in —iM?^^ and — , there are three logarithmic divergences hid- 

den in the UVDP parameter space, and they reproduce the corresponding sub divergences in 
the subdiagrams of Fig. 13(b). This feature was already anticipated by the electric circuits 
analogy of Feynman diagrams discussed in section III. However, when adopting a different 
set of UVDP parameters u, v, w transformed from the ones v%, t>2, t>3, the divergence regions 
in the parameter space are also changed correspondingly. 

Region (IV): —p 2 ^> m 2 . In this region, we obtain the first order correction to the wave 
function renormalization in the 4 theory. Clearly, there is no harmful divergence in this 
region as all the integrals of UVDP parameters are convergent. When — p 2 ^> m 2 , we can 
ignore all the terms proportional to m 2 in Ai 2 and the integral can be simplified to: 



(2)io 9 6(16tt 2 ) 2V y j J q (1 + w) 3 (1 + v) 3 y l + w 1 + 

duu ri M 2 

[ ln 5"^ 7 ~o\ - 7. + 1] 



[u + ( 1+M ,)( 1+1 ,)] 3 (l+u)[u(l+w)(l+v)+l] ( P 2 

i\ 2 2r l., M 2 

-p 2 [-(\n^-- lw + l) 



6(16tt 2 ) 2 ^ l 2 v -p 2 

+ Y^(-81 - 2 </>«(-) - 2 ^) + 2 ^ (1) (jj) + 2 ^«(~))], (114) 

where ip^(z) = ^lnT(z) is the polygamma function of order 1. 

Adding up all the contributions Eqs. (11051) . (11 10p (|TT3|) . fl 1 14[) . we arrive at the final result 
for the divergent contributions of the sunrise diagram Fig. 13(b) 

- iM % ] " {[3(ln 3 - ^ + 1]M2 - 3m2(ln ^ - 7 - )2 - 4m2(ln S - ^ 

+ ^ 2 ( ln ^|-7.)}- (115) 
2 — p z 

The counterterm diagram for Fig. 13(b) is shown in Fig. 13(61,2) and the result is simply 
given by: 



^ 2 V x ' J {2n) i k 2 -m 2 

i\ 2 M 2 M 2 

-2Ti^ (ln /- 7 - )[M2 - m2(ln #-^ + 1)] - (116) 

Note that we have used the vertex counterterm insertion of s and it-channels, so there is a 
factor of 2 in the above calculation. By summing up Figs. 13 (b) and (b\), we obtain: 

■ J i / r2(6)+(6i)+(6 2 ) ^A 2 fAf A/ r2 2\ 1 2/i M c \ , 1 2n M c \i 

~ %M {2) = (Io^^6 ( c _/i )_ 6 m ( 12 P ^^fl-^ 

+ lk» 2 - m 2 ) - im 2 (ln - \m 2 In ^ + ^ 2 In ^]} + - 
6 2 rrr 3 12 — ;r 

(117) 



As it is expected, the potentially harmful divergences m 2 (ln — 7^) In ^ cancel exactly. 
By considering the following overall counterterms for diagram (b): 



we get the final contribution to the two-loop self-energy as shown in Fig. 13(b): 

- <Mg» = jrf^lk^ 2 - m 2 ) - ^ 2 (ln 4) 2 - h 2 ln 4 + lV ln -4) + - (H9) 
1 ; (lovr^ 2 m z 3 m z 12 — 

We are now in the position to put all the results from diagrams Fig. 13 (a) and Fig. 13(b) 
together, and obtain the total contributions to the two-loop self-energy, 

iX2 A, 2 2m V 2 , 1/ 2 2n 



4 12 12 — 



Considering the massless limit m 2 — > and ignoring the quadratic contribution /x 2 — >■ 0, one 
arrives at 

" ^(2)K = 10/ !f ^ P 2 In A^, (121) 



which agrees with the one obtained by using the standard dimensional regularization 
method (see page 345 of the book[16|). 



C. Vertex Contribution at Two Loop 

The two-loop vertex contribution for the s-channel in 4 theory involves four groups of 
diagrams as shown in Fig. (I15p . It is expected that all harmful divergences cancel separately 
within each group. As the groups of diagrams (d) and (e) are related by a simple interchange 
of initial and final momenta, they should give the same results and it only needs to calculate 
either group and then to multiply by a factor of 2. 

First, let us calculate the simplest two-loop diagram (f) and its counterterm diagram (fl) 
in Fig. 14. Since the only UV divergence in (f) can be completely canceled by that in (fl), 
the group (f) does not require two-loop overall counterterm. 



f d 4 k 2 i f d% i 3 



^ " 2 J (2vr) 4 k 2 - m 2 J (2tt) 4 {k 2 - m 2 ) 2 ^ + p) 2 - m 2 } 

-> M 2 - m 2 In - 7^ + 1 , \ \ - \ dx 1 - x) 

2 16vr 2L c v m 2 ' 7J r(2)r(l) J v ; 

d*h 1 

(2tt) 4 {[(1 - x){k 2 - m 2 ) + x[{h + p) 2 - m 2 ]} 3 

tX 3 r , ,o n „ M 2 _ f 1 , 1 - x 



2(16vr 



2^2 



M 2 /" 1 — x 

[M| - m *(ln -| - 7 „ + 1)] | «fa m ,_ g(1 _^ . (122) 




FIG. 15: 



IA (2) -l[pO z m2 ){ IA) / „,2Wfc, ■ ^2 _ m 2 



(27r) 4 (A;2-m 2 ) 2 p 1 +p) 2 -m 2 
iA 3 r ,. - „ s o„ M 2 „ Z" 1 , 1 - x 



2(16vr 



2A2 



M 2 C 1 1 — T 

[(M 6 2 -/i 2 )-m 2 (ln^- 7oJ )] / dx— — , (123) 

J m — x{l — x)p z 



-'< )+(/1) = -^W[(^-- 2 )-- 2ln 4] f 1 dx ~ — ~7T~ — \~2 " ( 124 ) 
1 ; 2(167r / ) z m 2 Jo m — x (l — X )P 

The computation of diagram (c) and its counterterm diagrams (cl) and (c2) in Fig. 14 
is also straightforward, as these diagrams can be factored into a product of two one-loop 
integrals. The result is: 

-zAg = \(-i\) 3r lH ' 



W ~ 4 V ' 7 (2vr) 4 P-m 2 (A; + p) 2 -m 2j 
zA 3 r /■ d A k rl 



dx- 



4 7 (2tt) 4 J {(I - x)(k 2 - m 2 ) + x[(k + p) 2 - m 2 }}' 
*A 3 r f 1 . „ M 2 



dx(ln^ _- 7w )]^ (125) 



4(167r 2 ) 2 Jo m 2 — x(l — x)p 



and 



-aS )+(c2) =2.i(-^)HA) 1 <H 



( 2 ) ' 2 V AA V (27r) 4 A; 2 - m 2 (k + p) 2 - m 2 

iX 3 M 2 f 1 M 2 

= 2(m^ {ln 7# " 7 H dx(la m'-xd-xw ~ 7jL (126) 

where in the calculation of counterterm diagrams, the factor 2 in the first line accounts for 
the two equal diagrams (cl) and (c2). By adding up the above results, we have: 

. . (c)+(cl)+(c2) «A 3 f /i 2 2 M 2 2 

By considering the two- loop overall counterterm for diagram (c): 

" 4 (16vr 2 ) 2lln /i 2 7wJ ' ^ 
we obtain the renormalized contribution from the group (c) of diagrams: 

• \ 3 2 

-iA®=-—±——[ dx\n — £ — ] 2 . (128) 

(2)i? 4(16tt 2 ) 2 7 m 2 -x(l-x)p 2j v 7 

For diagram (d) in Fig. (fT5|) . due to the complicated dependence on external momenta, 
we shall focus on a simplified situation where only the s-channel contributes. Also, we only 
keep the leading divergent contributions, namely the log • log and log terms, because Ai 2 
can be simplified in such a subdivergence region. To see the momentum dependence in Fig. 
(|T5l) . we write explicitly the expression for diagram (d): 

% ( 2 ) ~ 2 [ % ' J (2tt) 4 J (2tt) 4 k\ - m 2 (k x + p) 2 - m 2 k 2 - m 2 (k 2 + k x + p 3 ) 2 - m 2 

iX 3 f d 4 k x f d A k 2 f 1 , 1 

dx 



2 J (2vr) 4 y (2vr) 4 y [x(h + p) 2 + (1 - x)k\ - m 2 } 2 
1 



\k\ — m 2 ] [{k 2 + k\ + J93) 2 — m 2 ] 
iX 3 f 1 . f d 4 k x f d A k 2 



dx 



2 j J (2vr) 4 J (2vr) 4 [k 2 + x(l - x)p 2 - m 2 ] 2 

' (129) 



[A; 2 - m, 2 ] [(/c 2 + k x + p 3 - sp) 5 



m 



21 



Note that in the process to obtain the second and third equalities, we have transformed the 
original integration into the general a/?7 integral with a = 7 = l,/3 = 2 and m\ = m 2 = 
m 2 , m 2 . = m 2 — x(l — x)p 2 by adopting the usual Feynman parametrization and making the 



translation of variable k% — > k% — xp. Thus, the general formulae Eq. (1751) gives, 

3 3 



\3 rl poo 6 6 i 



3=1 3=1 



.i 



(1 + v 2 ){3 + Vl +v 2 + v 3 ) 2 J (2tt) 4 [I 2 - M 2 } 2 

3 3 



/*1 poo J I 



i\ 3 

— > 



.1 

M 2 



with 



(l + v 2 )(3 + v 1 +v 2 + v 3 ) 2X M 2 



M 2 = m 2 _ X ( l - X )P 2 (P3 - X Vf 



(ln-JL- 7w ), (130) 



1 + V 2 3 + f i + f 2 + 1>3 

In order to carry out the above integral, we transform the UVDP parameters Vi,v 2 ,v 3 to 
the new set u,v,w as shown in Eq. (llOOp . so that the form of — «A/5 is changed to: 



( 2 ) 2 ■ (i6tt 2 ) 2 ,/ 7 (l + w) 2 (l + i;) 2 v 1 + w l + v J 

°° u M 2 . , 

r J; I i p^T^-'U ( 131 ) 



where .M is given by : 



.M 2 = m 2 -^*(1-*) P 2 " iP3 - Xpf 



u + 1 K ' u+1 u(l + w)(l + v) + 1' 

Notice that in Eq. (11311) only the integration over u is logarithmically divergent, while the 
ones over w,v convergent, so we can make the following approximation for M 2 in the limit 

u — > oo 

M 2 ~ m 2 - x(l - x)p 2 , (132) 
and the integration can be performed as follows 



( 2 ) 2-(167r 2 ) 2 7 y (1 + w) 2 (1 +v) 2 v 1 + w 1 + v 

°° « r M 2 

r i 1 19 2 /i \ 2 

[ M + (!+„)(!+,) ] m - < l ~ X )P 

i\ 3 f 1 M 2 M 2 

-+ -2T(T^T Jo dx{ln it~^ + mn m 2 -x(l-x) P 2 -^ (133) 

where the mass scale q 2 is taken to be m 2 — x(l — x)p 2 as the only scale in the limit u — > oo 
is Ai 2 ~ m 2 , = m 2 — x(l — x)p 2 . Thus we can write down the regularized expression for the 
diagram (d) as: 



{d) i\ 3 f 1 M 2 C 2 

'( 2 ) " ~2-(16tt 2 ) 2 J [{ m 2 - x(l - x)p 2 

^ rn'-S-W -^- (134) 



The counterterm diagram (d\) in Fig. 14 can also be simply calculated: 

_ iA w _ h^xM-t^) f i i 

1 W ~ 2 l A A J 7 (2vr)4P-m2(A; + p) 2 -m2 

^ 3 /"^ n M c u. Ml 

= ona 2 \2 / «^ m — ~7« m — 71 r-2-7uJ, (135) 

2(167r 2 ) 2 Jq // 2 m 2 — x(l — x)p l 

where 5^ accounts for the t— and u— channels in the subdiagram. By combining the above 
two equations Eq. (11341) and (1135)) . we have 

- * A S +(dl) * - o7^5 f 2 n , 2 ) 2 + ( ln 2 5 ri - 7 W )]- 

v J 2(lo7r z ) z Jo rrr — x(l — x)p z m z — x(l — x)p z 

Write the two-loop overall counterterm for the diagram (d) as 

• x (2)(d) «A 3 n \ Moci s 

-2(i6^ (ln ^- 7w) ' (136) 

then the two-loop vertex contribution of the diagram (d) is, 

((f) ^ 3 Z" 1 A* 2 2 V 2 

Obviously, the diagram (e) and its counterterm diagram (el) in Fig. 14 gives the same result 
as diagram (d) and (dl). 

By summing up the renormalized results of the diagrams (c), (d), (e) and (f) in Fig. (|15|) . 
we finally obtain the s-channel two-loop correction for the four-point function: 

-ikS'L = -J^{\[Cdx\jy. 

' 4 Jo 



m (16tt 2 ) 21 47o m 2 -x(l-x)p 2i 

+ C dx[{hx — £ — ) 2 + In /r> 



o 



m 2 — x(l — x)p 2 m 2 — x(l — x)p 2 

+ i[(/i 2 -m 2 )-m 2 ln4] / dx 2 1 ~ - J , (138) 
2 m z J m z — x{l — x)p z 

where — p 2 = —{pi + P2) 2 = s for s-channel. 

For completeness, we also have to consider t- and w-channel contributions, which will 
give us the similar expressions except for the substitution of p 2 : —p 2 = — (pi — P4) 2 = t 
for t— channel and — p 2 = — (pi — P3) 2 = u for u-channel. In Fig.( lT5]) . we only present 
counterpart of diagram (c) for t- and u-channels, and other diagrams in groups (c), (d), (e) 
and (f) could be obtained by the same change of external legs. 

As a consistency check for the above results, let us consider the situation with massless 
limit m — > and s — > 00 but keeping t fixed. From the identity s + t + u = 0, which indicates 
u ~ — s — > —00, one needs to consider u-channel as well. In this case, the four-point vertex 
correction is found to be: 

-zASt (u) ~ --^-^(In^f + ^ln^ 2 ! 



W " (16vr 2 ) 2L 4 v s> A' 



u 



5 i\ 3 fi 



2 



- 2 (i 6 ^< ln T> (139) 

where in the last equality we only keep the leading external momentum s = —p 2 ~ u 
dependence when the external momentum s is large. Such a result agrees with the one given 
in the book[l6|(on page 345). 




FIG. 16: 

D. Two Loop (3 Functions and Anomalous Mass Dimension in </> 4 Theory 

Having calculated the divergence behavior of all the two-loop diagrams in the </> 4 theory, 
we are now ready to obtain two-loop f3 functions. From Eqs. (1961) and (11181) . we can extract 
the two-loop mass and wave function counterterms respectively: 



-«3 = ^^ M "-^f -7.)-i(M c 2 -/, 2 ) 



1 ,2^ M c „ ^2 , 1 2. M c 



- m'Qn -§■ - lu ) A + -m 2 (ln -f - 7w )], (140) 

And the two- loop vertex counterterm can be extracted from Eqs. (11271) and (11361) : 

i\ 3 . 3 n M 2 , 2 0/1 M c 2 ,_ n , 

= (ifaf '"I 1 "i# -7 "' + ' (142) 

where the factor 3 in the first line accounts for the s, t, u-channels respectively. 

Recall the relation between renormalized coupling constant A and the bare one Ao is: 

A = A Z| - 5\ = A (l + 5 Z ) 2 - 5 X 

« A (l + 26 z ) - S x , (143) 

where 5z and 5\ are the function of the bare coupling A , which is independent of scale \i. 
In the perturbative calculation of A at two-loop level, we have, 

AS3A °- 2 'l2(fe< h f-'»> 

3A ° ,inM_ 7 „ ) _ 7T ^ [ 3 (ln M_ 7 „ )2 _ 3(ln M_ 7 „ )] 



2-(16vr 2 ) v /z 2 m ' (16vr 2 ) 2L 4 v /i 2 v // 

A ° " Mi% (ln / - ^ - (i6# [ 4 (ln / - ^ - ¥ (ln / - ^ )] ' (144) 



Thus according to the definition of /3-function which is supposed to sum up all the leading 
logarithmic terms (ignoring the logarithmic-squared term), we arrive at the /3-function for 
the renormalized coupling constant A as: 



dX 

Pa = 

3A 2 



16vr 2 

3A 2 

16tt 2 



2X1 



17 



(16tt 2 ) 2 6 
17 A 3 

2^2 : 



3 (16vr 



(145) 



where the bare constant Aq has been replaced in the last line by its renormalized one, leading 



to the standard result (3\. [18H20 



Similarly, we can evaluate the anomalous mass dimension at two-loop level. From the 
definition of the renormalized mass: 



m 



Z^m\ - 5 m 2 = ml + m 2 5 z - 5 n 



(146) 



we have the following approximate relation for the renormalized m 2 given in terms of bare 
mass m 2 , and the bare coupling constant Aq at two-loop level: 



m 



m 2 + 



A 



2(16tt 2 



M 2 

■[(M c 2 -/i 2 )-m^(ln-f - 7w )] 



\2 
A 



1 



M. 



7.)-^(M c 2 -^ 2 ) 
Ml 



-m (ln— Y-n/w) + — m (ln— ^ - 7w )] 



(147) 



which is different from the result obtained by using the dimensional regularization approach 
due to the appearance of the quadratic terms. The anomalous mass dimension where we 
sum up all the leading quadratic and logarithmic terms (i.e., not considering the logarithmic- 
squared term and quadratic-logarithmic cross term) is given by: 
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;i67r 2 )m 2V 2 
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(148) 



where we have replaced in the third line the bare mass and coupling constant with the 
renormalized ones. 

Note that the resulting 7^,2 is different from that obtained in ref. [2l[ by using the dimen- 
sional regularization approach with the MS subtraction scheme. The difference occurs in 
both the power-law and the logarithmic running terms. For the power-law running terms 



with the form n 2 /m 2 in 7^2, it reflects the fact that the LORE method maintains the origi- 
nal quadratic divergence. For the logarithmic terms, the difference can be caused from the 
well-known fact that the two-loop anomalous mass dimension in 4 theory is in general sub- 
traction scheme dependent. This may be seen from the rescaling /x 2 — > e°°/x 2 , the resulting 
leading logarithmic term at two loop level is changed by an additional contribution from the 
logarithmic-squared term, thus the corresponding 7^2 for the logarithmic running is changed 
to be 



As a consequence, both the /i 2 -independent term and the quadratic /^-dependent terms also 
changed correspondingly. Similarly, when shifting the scale /1 2 — > p? = fi 2 — aom 2 , the 
leading logarithmic term also receives an extra contribution from the quadratic-logarithmic 
cross term, and the resulting 7^2 for the logarithmic running in terms of the new subtraction 
energy scale fi 2 is modified to be 



However, the quadratic-logarithmic cross term is now given in terms of two energy scales /1 2 
and fi 2 rather than a single one, i.e., (M 2 — /i 2 )(ln M 2 /fi 2 — 7^). From the above illustration, 
it is seen that either the rescaling or the shifting of the subtracted energy scale /z 2 will change 
the initial correlative form (M 2 — /i 2 ) and lnM 2 //x 2 . Therefore, when the quadratic terms 
are kept by using the LORE method, the arbitrariness caused by the subtraction scheme 
for the scalar mass renormalization at high loop order may be eliminated by requiring to 
maintain the correlative form (M 2 — fi 2 ) and lnM 2 //z 2 with a single subtracted energy scale. 

IX. GENERAL PROCEDURE OF LORE METHOD 

With the explicit calculations of two-loop Feynman diagrams in the 4 theory given 
above, it is useful to summariz the general procedure in applying the LORE method to 
multi-loop calculations. It is expected that the same procedure is applicable to higher-order 
calculations with similar features when merging the LORE method with the Bjorken-Drell's 
analogy between the Feynman diagrams and electrical circuits, though we have only shown 
it in the two-loop calculations. The procedure may be wtated in the following steps: 

(i) Write down the corresponding Feynman integrals by using the Feynman rules of the 
theory for any given Feynman diagrams. 

(ii) Combine the denominators by using Feynman parameters to evaluate the two-loop 
integrals into the sum of the a /?7 integrals of scalar- type and tensor-type. The use of 
the usual Feynman parametrization in this step needs to be distinguished from the UVDP 
parametrization adopted for the a/?7 integrals. The latter may contain the UV divergences, 
while the former is in general irrelevant to the UV divergences but it can contain infrared (IR) 
divergences. From this point of view, making distinction of Feynman parameters from UVDP 
parameters enables us to separate IR divergences from UV divergences in two parameter 
spaces. 

(iii) By applying the general formulae Eqs. (I75H77P for the ILIs of two-loop a/3j integrals to 
the resulting aftj integrals coming from a given Feynman integrals, we can straightforwardly 
read off the final results for those integrals. Alternatively, one may also adopt a practically 
useful procedure by completing the squares of the factors in the denominator and evaluate 





the a/?7 integrals into the 2-fold ILIs as proposed in ref. [3|], which shows that for each internal 
loop momentum, one can always transform the integrals into the 1-fold ILIs with respect to 
it, and then integrate out the ILIs by means of the LORE method. The two procedures are 
actually equivalent. For tensor-type integrals, we need to apply the consistency conditions 
Eq. (jT3l) to transform them into the corresponding scalar- type ones first. 

(iv) When the integrals involve overlapping divergences, the above procedure will trans- 
form the divergences appearing in the subdiagrams into the ones in the UVDP parameter 
space. In order to identify those divergences, it is helpful to use the advantage of the 
Bjorken-Drell's analogy between the Feynman diagrams and electric circuits. To extract the 
UV divergence behavior, it is useful to explore the possible divergence regions in the UVDP 
parameter space. Then apply the prescription described in Eq. (1321) through introducing a 
mass scale q 2 to transform the integrals into the momentum-like ones, so that we can directly 
apply the LORE method. The scale q 2 Q is in general taken to be the renomalization scale 
or some intrinsic scales in the original Feynman integrals, such as the masses of particles 
and/or the external momenta. The explicit form of q 2 Q should be fixed by certain criteria, 
such as the typical scale in the divergent regions of parameter space, so that the harmful 
divergences cancel exactly. 

It is interesting to notice that in the LORE method only the overall divergence of the 
overlapping Feynman diagrams is expressed in the momentum integration and the resulting 
functions yo(x) and yi [x) can depend on the mass factor M 2 through which a dependence 
on kinematic invariants comes in. In contrast, all other divergences arising from the sub- 
diagrams are actually given in terms of the UVDP parameters and the resulting function 
y$ in the logarithmic divergence will be independent of any kinematic invariants. This can 
be seen from the general scalar- type integral Eq. (1611) which can arise from n-loop Feynman 
diagrams. The more detailed evaluations are carried out for the so-called a/?7 integral of 
two loop diagrams, which can explicitly be seen in Eqs. (l77tl84|) . As a consequence, the regu- 
larized divergent quantity of subdiagrams involves only a kinematic-independent scale fi via 
a polynomial of fi 2 /M 2 with /j, 2 /M 2 — > at M c — > oo, thus the function yo(x) arising from 
the subdiagrams is no longer a complicated function of kinematic invariants. 

X. CONCLUSIONS AND REMARKS 

We have explicitly shown how the loop regularization (LORE) method can be consistently 
applied to two loop calculations of Feynman diagrams, appropriately treating the overlap- 
ping divergences. The key concept of the LORE method 0, H| is the introduction of the 
irreducible loop integrals(ILIs), which are generally evaluated from the Feynman diagrams 
by using the Feynman parametrization and the ultraviolet-divergence-preserving(UVDP) 
parametrization. We have demonstrated in this paper how the evaluation of ILIs and UVDP 
parametrization naturally merges with the Bjorken-Drell's analogy between Feynman dia- 
grams and electric circuits. In particular, the UVDP parameters can be regarded as the 
conductance or resistance in the electric circuit analogy, and the sets of conditions required 
for evaluating the ILIs and the momentum conservations have been found to associate with 
the conservations of electric voltages in each loop and the conservations of electric currents 
at each vertex respectively. As a consequence, the divergences in Feynman diagrams corre- 
spond to infinite conductances or zero resistances in electric circuits, and the LORE method 
merging with the Bjorken-Drell's analogy has the advantage in analyzing the complicated 
overlapping divergence structure of Feynman diagrams. Therefore, the Bjorken-Drell's cir- 
cuit analogy allows us to clarify the origin of UV divergences in the UVDP parameter space 



and identify the correspondence of the divergences between sub diagrams and UVDP param- 
eters. From the explicit calculations of the case with a = /3 = 7 = lin the general aj3^y 
integral, the divergences arising from the subintegrals manifest themselves in the integration 
over the corresponding asymptotic regions of the UVDP parameter space. The calculations 
of the corresponding counterterm diagrams confirm our intuitive picture that all the harmful 
divergences cancel exactly in the final result. Although the procedures and calculations in 
the LORE method are not as concise as the ones in the dimensional regularization, the over- 
lapping divergent structure and behavior as well as its treatment become more physically 
clear in the LORE method. 

As an interesting application, we have taken the massive scalar 4 theory as an example 
and performed the detailed calculation of two loop contributions by applying the general 
formalism of the LORE method. By explicitly computing the two- and four-point functions 
at two-loop level and carefully using the advantage of Bjorken-Drell's circuit analogy, all 
the harmful divergences cancel exactly and the resulting two loop corrections agree with the 
standard results for the logarithmic corrections. The power-law running of mass is explicitly 
given at two loop level. 

In this paper, we have only carried out two-loop calculations and explicitly demonstrated 
the consistency of the LORE method at two-loop level, but it can be shown that the general 
procedure of the LORE method shown in Eqs. (I57H6T]) is applicable to even higher-loop 
calculations by taking advantage of Bjorken-Drell's circuit analogy. Furthermore, we only 
considered the scalar-type two-loop integrals. However, as shown in [4], in order to ensure the 
gauge invariance, it is necessary to keep the consistency conditions Eq. (fl3l) which correctly 
transform the tensor-type ILIs into the scalar-type ones. We shall demonstrate how these 
consistency conditions in two-loop or even higher-loop order by an explicit calculation j22|. 
although it has already been demonstrated in a general way inj^]. We would like to point 
out that the advantage of merging the LORE method with Bjorken-Drell's circuit analogy 
enables us to figure out a more general and rigorous proof for the validity of the LORE 
method to all orders in the perturbation theory [231] . 
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Appendix A: Useful Formulae in UV Divergence Preserving (UVDP) Parametriza- 
tion 

The introduction of UVDP parameters is to combine the various denominators propagat- 
ing factors, whose utility is similar to Feynman parameters. The motivation to introduce 
a new UVDP parametrization method is to transform a divergent integral in the UVDP 
parameter space into a ILI-like divergent one, the object regularized by the LORE method. 
The simplest case is to combine only two factors in the denominator by using the identity: 

1 f°° du dv 1 1 1 

ab~J (i+n) 2 (i+^) 2 ( "r+^r+v^ + ^p- 



If one of the factors have more than one power, we can differentiate with respect to A or B 
to get, 



1 /' x (In ilr ^ 1 1 N JT+^y 



AB- J (1 + «)3(1 + V )2 ^ 1 + u 1 + V > + _J_]n+r 

More general identity for more than two factors is: 



1 



r fr d Vi ^ 1 _ , (n-1)! 

io fia+^) 2 j ELAi^ 



jM 2 • • • A Ti 
Even more general form can be derived: 

i r°A „,a i _ nr=i 



Alternatively, we may also take another more useful form for the case of two factors by 
just integrating out one of the parameters u and v by using the delta function, which has 
been adopted in [jjj]: 

i r° , nu n - 1 

du [A , „, D1 ^ 1 ; (A2) 



but this form cannot be generalized to the more general case easily. 



/,From the general identity Eq.flATJ), we notice that the relation of the UVDP parameters 
Vi to Feynman parameters Xi is: 

Xl = — |— . (A3) 

This identification allows us to transform a divergent integral with Feynman parameters into 
the one with UVDP parameters, which can be further transformed into a ILI-like integral by 
introducing a free mass scale and being regularized in the framework of the LORE method. 
Such a trick is discussed in Eq.( 
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